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Abstract 

The longwavelength excitations of a quark-gluon plasma at high temperature 
can be described as collective oscillations of gauge and fermionic average fields. We 
show that, at leading order in the coupling strength, the Dyson-Schwinger equations 
for the iV-point functions reduce to a set of coupled equations for these average fields 
and their induced sources which involve only 2-point functions. The equations for 
the 2-point functions describe the dynamics of the hard plasma particles in the 
presence of soft background fields. They may be given the form of simple kinetic 
equations. Both the wavelength and the amplitude of the collective modes are 
controlled by the coupling strength, and we show that it is important to take this 
properly into account in order to obtain consistent equations of motion which are 
covariant under gauge transformations. By solving the kinetic equations for the 2- 
point functions with well defined (i.e. retarded or advanced) boundary conditions, 
we obtain the induced currents in closed forms. These act as generating functionals 
for the one-particle irreducible amplitudes with soft external lines, and yield in 
particular all the so called "hard thermal loops" identified in the diagrammatic 
analysis. 
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1 Introduction 



The study of the long wavelength excitations of the quark-gluon plasma is interesting 
in several respects. In particular, it provides clues to basic features of the plasma off- 
equilibrium dynamics, and it is a first step towards a full description of its transport prop- 
erties. Interesting progress in this area has been made recently, and consistent schemes 
to perform calculations of amplitudes in hot gauge theories have begun to emerge [1-6], 
giving rise to various applications [7-11]. 

We consider in this paper a quark-gluon plasma in thermal equilibrum at high 
temperature T well beyond the critical temperature for the deconfining phase transition, 
and neglect quark masses. Then, the only energy scale in the problem is T. We also assume 
a regime where the effective coupling constant g is small. In the plasma, quarks and gluons 
have typical energies and momenta of order T, and will be referred to as "hard" particles; 
their number density is proportional to T 3 . Thus, the typical interparticle distance is of 
the order 1/T, that is, of the order of the thermal wavelength At- This means, for instance, 
that the Pauli principle cannot be ignored for the hard particles. However, when coupled 
to weak and slowly varying disturbances, the plasma may develop a collective behaviour 
on a typical length/time scale of order 1/gT. A familiar example of such a collective 
behaviour is Debye's screening. Collectivity arises because any motion taking place over 
a distance scale 1/gT may involve coherently a large number of hard particles. We shall 
see later that such collective motions can be treated as classical degrees of freedom. 

The existence of such possible excitations manifests itself in specific difficulties en- 
countered in perturbative calculations. In particular, diagrammatic analyses have revealed 
that amplitudes with "soft" external lines, i.e. with all external momenta of order gT, 
may receive contributions from an infinite number of loop diagrams which are as large as 
the tree amplitudes [O, [TJ|, ELR EL EL These contributions arise from one- loop integrals 



in which the loop momentum is hard, and for this reason they have been named "hard 
thermal loops" (HTL). As they are of the same order of magnitude as the tree amplitudes, 
they have to be resummed in consistent higher order calculations JTJ, |2|] . 

The need for resummation in a perturbative expansion usually points to simple un- 
derlying physics, as shown, for instance, by the classic example of the electron gas. There, 
the resummation of ring diagrams performed by Brueckner and Gell-Mann |20| reflects 



the simple physical phenomenon of Debye's screening. In order to illuminate the simple 



physics underlying the resummation needed in QCD, we have initiated JT6], |17]] a new ap- 
proach to the problem, by working directly at the level of the equation of motion rather 
than on the Feynman diagrams. The crucial remark is that the collective phenomena that 
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we are trying to describe involve mean fields which represent the average motion of hard 
particles over distances large compared to their mean separation. Such mean fields may 
be induced by external disturbances, but they are also the natural degrees of freedom 
describing the soft normal modes of the quark-gluon plasma. Our main objective will 
be therefore to obtain consistent equations of motions for these mean fields. In order to 
allow for collective excitations with arbitrary quantum numbers, we shall consider both 
fermionic and bosonic average fields. The physical picture which will emerge from our 
analysis is that the dominant interactions in the plasma are those between the mean fields 
and the hard particles. Direct interactions among hard particles can be neglected, but 
since hard particles can contribute coherently, they strongly affect the dynamics of the 
mean fields. Our final equations will take the form of a coupled set of mean field equa- 
tions and kinetic equations describing the motion of hard particles in the presence of the 
mean fields. These equations generalize the familiar set of Vlasov and field equations of 
ordinary electromagnetic plasmas |L8], [0| . 

When dealing with this type of problems one is usually facing three kinds of ap- 
proximations, involving, respectively, the strength of the coupling, the amplitude of the 
oscillation, and the wavelength or period of the modes that one is studying. In the present 
case, there is no other scale than the temperature, and it turns out that all three approx- 
imations are controlled by the single small parameter, g. They should be viewed, then, 
as three facets of the same approximation, namely an expansion in powers of g. Thus, 
the collective modes with wavelength ~ ^/gT, have, as we shall see, amplitudes limited 
by i[)ip ^ gT 3 and ^ gT 2 , respectively. In order to preserve gauge covariance, it is 
essential that these three features of the approximation scheme be implemented together 
consistently. The naive loop expansion fails precisely because it takes care only of the 
first facet, namely the expansion in powers of the interaction vertices. 

As alluded to before, one outcome of the present analysis is a set of gauge covariant 
kinetic equations for hard quarks and gluons moving in soft gauge and fermionic fields. 
We thus make contact with previous works on the kinetic theory of non abelian plasmas 
[21-25]. However, the special role of the scale gT, as a unifying parameter controlling the 
various approximations being done in deriving kinetic equations, has never been properly 
recognized. Because they include terms which are found to be of different orders in our 
expansion scheme, most kinetic equations proposed so far are not only more complicated, 
they are also inconsistent. Besides, although quite some effort was spent in using gauge 
invariant Wigner functions, specific issues related to gauge invariance in a quantum theory, 
such as gauge fixing independence, or the elimination of ghost particles, have never been 
properly addressed. The resulting inconsistencies can be seen in calculating one-particle 
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irreducible amplitudes which, in leading order, can be deduced from the kinetic equations. 
As well known, these amplitudes satisfy Ward identities imposed by gauge symmetry. This 
puts constraints on the structure of the non-linear terms which enter the kinetic equations, 
and such constraints do not appear to be satisfied by the transport equations proposed 
so far. 

Our strategy for deriving the relevant equations of motion is similar to that already 
used in connection with the abelian case However, the structure of the non abelian 
theory is much richer, and the present work contains many new technical developments. 
It is organized as follows. In Sec. 2, after summarizing the conventions and the notations 
used in this paper, we write the equations of motion relating the average fields to the 
off-equilibrium 2-point functions. The connected parts of these 2-point functions play 
the role of induced sources for the average fields. From them, one-particle irreducible 
amplitudes for the soft fields can be calculated. 

The equations for the average fields constitute the first level in the hierarchy of 
Dyson- Schwinger equations for the iV-point functions. In Sec. 3, we describe all the ap- 
proximations which allow us to calculate the 2-point functions. In leading order, these 
functions concern only hard particles, for which radiative corrections can be neglected. 
This implies that the hierarchy truncates at the level of 2-point functions. We analyti- 
cally continue the 2-point functions to real time, and we implement the condition that 
the average fields are slowly varying with the help of Wigner functions and gradient ex- 
pansions. We also take into account the fact that the strength of the mean fields have 
to be constrained so as to obtain consistent equations. We emphasize that this does not 
result in a trivial linearization of the equations since, for example, gauge fields involved 
in covariant derivatives are to be kept to all orders. At the end of Sec. 3, we have a set of 
gauge covariant equations for the induced pieces of the 2-point functions. These equations 
incorporate consistently all leading order effects in g. 

In Sec. 4, the equations for the 2-point functions are transformed into equations for 
the induced sources. These describe the polarizability of the plasma at leading order 
in g. They involve only the physical hard particles from the thermal bath, i.e. the 
quarks, the antiquarks and the transverse gluons. The equations are covariant under 
gauge transformations of the mean fields. They reduce to simple kinetic equations for 
density matrices which, in general, are non diagonal in color space, or which can mix 
boson and fermion degrees of freedom. These density matrices play the role of generalized 
classical distribution functions for the hard particles. The kinetic equations are solved for 
retarded boundary conditions. We get then explicit expressions for the induced sources 
in terms of the average fields. 
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The knowledge of these induced sources entirely determines the dynamics of the 
average fields, which is discussed in Sec. 5. There we obtain all one-particle irreducible 
amplitudes for soft fields by functional differentiation of the induced sources. We recover 
in this way all the "hard thermal loops", with the boundary conditions appropriate to 
the problem at hand, i.e. the causal propagation of the mean fields. The Ward identities 
satisfied by the HTL are seen to follow simply from the conservation laws obeyed by the 
induced currents. In cases where Landau damping is inoperative, a generating functional 
for the induced currents is obtained. It may be regarded as the effective action which 
describes the propagation of long wavelength fields in a polarizable non abelian plasma, 
and coincides with the generating functional for hard thermal loops previously derived in 



26| . [27], |28| . Properties of this generating functional are also discussed in |29j, |30 



The main conclusions are summarized in Sec. 6. Technical details related to the use 
of a general covariant gauge in our main discussion are gathered in Appendices A and B. 
In App.A we show that, in leading order, ghost particles do not undergo collective motion. 
This is in line with our conclusions concerning the semiclassical character of the leading 
order polarization phenomena. This same conclusion is supported by the independence of 
our results with respect to the choice of the gauge in the original quantum theory. This 
is proved in App.B. In App.C, we clarify the conditions under which the effective action 
proposed in Sec. 5 generates indeed the renormalized mean fields equations. 



2 Mean fields and induced sources 

In this section, after specifying some of the conventions and the notations used in this 
work (in subsection 2.1), we derive the first set of equations in the hierarchy of Dyson- 
Schwinger equations for iV-point functions, i.e., the equations of motion for the average 
gauge and fermionic fields (in subsection 2.2). These equations relate the mean fields 
to 2 and 3-point functions which are interpreted here as induced sources for the mean 
fields. Functional relations between these induced sources and the one-particle irreducible 
amplitudes are recalled in subsection 2.3. The relation between induced sources and 2- 
point functions is further analysed in the last subsection. Most of the material presented 
in this section is standard, and the equations written are exact; the approximations to be 
made will be presented in the next section. 
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2.1 Generalities 



We consider an ultrarelativistic plasma in, or close to, thermal equilibrum, at a temper- 
ature T = 1/(3. The equations of motion which describe its long wavelength excitations 
will be derived using the imaginary time formalism, and their real time versions will be 
obtained through an analytic continuation. Since the only quantities to be continued in 
the present work are at most 2-point functions, the procedure is non ambiguous. It allows 
us to avoid the intricacies of real time formalisms and to preserve, at least in the first stage 
of our approximation scheme, the familiar structure of perturbation theory for systems 
in equilibrium. In doing the analytic continuation, the relevant 2-point functions may be 
chosen to obey retarded boundary conditions, appropriate to the calculation of the causal 
response functions that one is eventually interested in. We use Minkovski metric and we 
note x = (xq, x) = (—it, x), 8q = id T and d 2 = —d 2 — V 2 . In the analytic continuation, xq 
becomes a complex variable (r acquires an imaginary part) and the metric, the gamma 
matrices, etc., are untouched. Our units are such that H — c — 1. We consider a SU(iV) 
gauge theory with Nf flavors of quarks . The color indices for the adjoint representation, 
a,b,... run from 1 to N 2 -l, while those for the fundamental representation, run 
from 1 to N. We work in four space-time dimensions, d — 4, but on occasion we shall 
keep d as an explicit parameter in our equations in order to exhibit the contribution of 
the various gauge degrees of freedom. Minkovski indices \i, v, ... run from to d — 1. 
The generators of the gauge group in different representations are taken to be Hermitian 
and traceless. They are denoted by t a and T a , respectively, for the fundamental and the 
adjoint representations, and are normalized such that 



Furthermore, (T a ) bc = -if abc , and t a t a = C f , where C f = (N 2 - 1)/(2N) is the Casimir 
of the fundamental representation and f abc are the structure constants of the group: 



We use, without distinction, upper and lower positions for the color indices. 

The thermal averages of fields operators in the presence of external sources can be 
obtained from the generating functional 






[t a , t b ] = if abc t c . 



(2.2) 




x exp-{s + (j,A) + (f},ij) + $,r j ) + (C,0 + (C,C)}, (2.3) 



5 



where S is the usual QCD action in a covariant gauge 

S = £ drSx + -(<?■ A a f + <AHW + C a ^(^C) a } , (2.4) 

and ip, i>, A £ and £ denote, respectively, the quark fields, the gauge fields and the Fadeev- 
Popov ghosts. They obey antiperiodic (for ip, ip) and periodic (for A, ( and () boundary 
conditions, with period (3. In the formula above, A is the gauge fixing parameter, and 
jp = with the usual notations for the covariant derivatives in the fundamental (D^) 

and in the adjoint representations {Dp): 

Dp^dp + igA^x), Dp^dp + igA^x). (2.5) 

Here A^ = A^t a and A^ = A^T a are the gauge field matrices in the two representations. 
The source terms in eq. (|2.3| ) are 

(j,A) = jyrd 3 x 3 ;(x)A^x), (2.6) 

and similarly for (fj,ip), (-0, 77) , (C, () and (£, C). 

The gauge field strength tensor is F^ v = [D^, D u ]/{ig) = F® u t a . For any operator 
of the form O(x) = O a (x)t a , we have 

[D„ 0{x)\ = d,0(x) + tg[A,(x), 0(x)}, (2.7) 

and similarly in the adjoint representation. 

2.2 Equations of motion for the mean fields 

In the absence of the external sources, the plasma is in thermal equilibrium and the expec- 
tation values of the fields vanish. For nonvanishing sources, the fields acquire nontrivial 
expectation values which obey equations of motion that are readily derived from eq. (|2.3|) : 

([D\F Ufl (x)}) + \d,d u {A v (x)) -gt a $(x) lfl t a i/;(x)) + gt a f abc {(d,t\x))C(x))= 3 ,(x), 

(2.8) 

tf{il>(x))-g(4(x)i/>(x)) = n(x), (2.9) 
- d 2 (C(x)) + gf abc (&* (A b ,(x)C(x)) ) = C a (x), (2.10) 
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- d 2 (C(x)) + gr hc (Al{x)d»C{x)) = C a (x), (2.11) 

together with the Hermitian conjugate equation for (ijj). The explicit form of the com- 
mutator term in eq. (|2.8|) is 

([D\ F v ,\ a ) = (g, u d 2 - dpdj) (A v a ) - gf abc T, pXu (d x A?) ) - g 2 f abc f cde (j^Af^ ) , 

(2.12) 

where 

^fiupx = ^9fiu9pX — 9fi P 9uX — 9p\9up- (2-13) 
Note the following identities which will be used later: 

k^T pV p\k = k v kp g V pk , 

^pupX ^uppX ^puXp r pXpuVi r jiupX ~\~ r^jpAi^ ^pXup 0. (2.14) 

We shall also use the following shorthand notation for the four gluon vertex: 

^tvXp = 2 { f eab f eCd (9p,x9u P — 9fj,p9ux) + f eac f edb {9np9xu — 9pu9xp) 

+ r d r h \9»v9px-9»xgp„)\ (2.15) 

Note that T ab ^ x an d r \ vp x ar e related by 

f^abcd Ap /\X rade rbce-p AP A x (O 1 P,\ 

1 pupX^c^d = ~J J 1 ppvX-n-c-n-d- l Z - iD J 

We separate the 2-point functions into connected parts, denoted by the symbol 
(■ • -} c , and disconnected parts. For example, (A(x)ip(x)) = (A(x)tp(x)) c + (A(x))(ip(x)) , 
and (ijj'y f j i t a "ip) = (ijj'-f tl t a ifj) c + (ip) / j fJj t a (ip). Then eq. ( [2.1 2|) can be written compactly as 
follows: 



{[D v ,F v »] a ) = [V\^] a - gr bc Tp pXv (Al(d x Ai)) c 

- 9 2 f abc f cde ( A\A%A v e ) c - g 2 T a ^ x ( Al) ( N> C A\ ) ^ , (2.17) 

where £> M = <9 M + ig(A fl (x)) and = [D^V v ]/{ig). 

After the separation of the connected parts in the 2 and 3-point functions, the mean 
field equations ([2.8|)-( f2~Tl| ) can be written as: 



[V\ T Vil {x)\ a + Xd,d u {Al(x)) - g^x))^^)) + gf abc (d,(C b (x))) (C(x)) 

= r^)+r; da (x), (2.18) 



(0 - g(A{x))) = -nix) + V md (x), (2.19) 

- d 2 (( a {x)) + gf abc d^ ((Afa)) (C{x))) = C a {x) + C inda {x), (2.20) 

- d 2 (C{x)) + gf abc {Al{x))d^{C{x)) = C a {x) + C mda (x). (2.21) 

In the left hand sides of these equations we have isolated the terms corresponding to the 
tree level equations of motion generated by the QCD action (|2.4j) . In their right hand 
sides, we have, next to the external sources, connected 2 and 3-point Green's functions 
which we refer as "induced sources", because they play a role similar to that of r], j M , or 
C, that is, they act as extra sources for the mean fields. Explicitly, they are 

f; da (x) = g(^t a ^) c - gf abc ((d^ b )C) c + gf abc T, pXu ( A v h (d x A^ ) c 

+g 2 f a X d x ( Aft ( A P C A\ ) + g 2 f abc f cde ( AlA%A\ ) , (2.22) 



,2 1 vilx-il / I i/» | \ ^ | g2 jabc jcde / j^l > j^d J^u 

c \ v ^ e I c 

rf nd {x)= gYt a {A a v {x)^{x)) c , (2.23) 
Cr M {x) = -gf abc (cP (Al(x)C(x)) } c , (2.24) 

C amd {x) = -gf abc (Afa^Cix)) c . (2.25) 

In equilibrium, when the external sources vanish, the induced sources are also vanishing. 
This follows from symmetry. In equilibrium, the expectation values involve thermal av- 
erages over color singlet states, and elementary group theory can then be used to prove 
that, in this case, all terms on the r.h.s. of eqs. (|2.22|) — (|2.25|) indeed vanish. 

At this point, it is perhaps useful to illustrate the structure of the equations of 
motion by means of diagrams. This is done in Figs.l and 2, which represent schematically 
the various contributions to the expectation values of the fields, as obtained from the 
equations of motion (|2.18|) and fl2.19|). 



2.3 Induced sources and one particle irreducible amplitudes 

It will be convenient later to express the induced sources as functionals of the mean 
fields. Following a standard procedure, this may be achieved in a systematic fashion by 
introducing the effective action, that is, the Legendre transform of InZ, with Z given 
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by eq. ( |2.3| ). The external sources are then considered as the first derivatives of the 
effective action with respect to the mean fields, which are treated as independent variables. 
Further derivations yield the one particle irreducible amplitudes. Thus, by differentiating 
eqs. (|2.18|) and ( [2.19| ) with respect to the mean fields, one gets, using the usual definitions 
of self energies, 

Si a (x) 8i inda (x) 

S^ v + Ut(x,y), nj^y)^, (2.26) 



5(A" b (y)) u ^ ^ ^ 5(Al{y)) 



,)//( ' r) ' i$) + g(A(x)))5(x-y) + X(x,y), ^ y ) = S J^M. ( 2 .27) 



5(^{y)) v r ™ v v " >"> v ' y/ 5(^{y)) 

In eq. ( ggg) 

is the free inverse propagator for gluons. A further differentiation gives the 3-point func- 
tions. The quark-gluon vertex is 

trf? + aTHx v z) = 52] " {X) = ^ (2 29) 

91 +9L « [X > V > Z) 6Mz))6Mv)) 6{A*{x))6Mz)y ^ 

that is, 

c2 -inda 



8Z(y,z) 5% n - [x) 



gr%z, y, z) = 77^^ = FTTT^TTTTbv ( 2 -30) 



Similarly, we get for the vertex correction to the three-gluons amplitude: 



S(Al{y))S(AUz)) S{Al(z)) 



We do not write here the expressions corresponding to amplitudes with ghost external 
lines since we shall not need them. 



2.4 Induced sources and two point functions 



By inspecting eqs. (gjgj)-(g25D, one sees that the induced sources are given by off- 
equilibrum 2-point functions. The only exception is the induced current j md which also 



contains the 3-point function ^ A^A^A" J . However, in the present work, we shall be 
interested only in the contributions to the induced sources which appear at leading order 
in g. Since the connected 3-point function contains at least two powers of g more than 
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the other terms, it can be safely ignored. This allows us to express the induced current 
jind ( |2.22 ) in terms of 2-point functions only. We write 

f; da (x) = + f G ,(x)+f g ,(x)+f n ,(x), (2.32) 

where the subscripts (f, G, g, Q) refer to which hard particles contribute to the induced 
current. Thus jf is the current associated with the polarization of the hard fermions 
induced by the average quark and gauge fields 

j^(x) ee ji^/^^c (2.33) 

Similarly, j G results from the polarization of the hard ghost particles induced by the gauge 
mean fields 

f Gfl (x) = -gf abc ((d,C b (x))C(x)) c . (2.34) 

The role of this current is to eliminate spurious polarization effects which, in a covariant 
gauge, are due to non-physical gluonic degrees of freedom. The induced current resulting 
from the polarization of the hard gluons by the quark and gauge mean fields is 

= 9f abc r w x» ( Aftx) (d x A?(x)) ) c . (2.35) 

The last term in eq. (|2.32|) is the "tadpole" contribution 

jW = 9 2 ^%{^(x)} (A» c {x)A\{x)) c . (2.36) 

It is convenient to express the induced sources in terms of time-ordered propagators, 
for which we shall derive equations of motion in the next section. Let us first introduce 
specific notations. Aside from the usual quark, gluon and ghost propagators, 

S(x,y) = (T4>(x)^(y)) c = ' ){r( ' r)} 



5r](y) 

Df v {x,y) ee {T A;(x)Al(y)) c =- 5 -^^-. 



G ab (x,y) ee (TC(x)C(y)) c = -^kp ( 2 - 37 ) 



we shall also need the following "abnormal" propagators 
K a v {x,y) ee (T^(x)A:(y)) c 

H a u (x,y) ee (TA a u (x)$(y)} c 

O ab (x,y) ee {TC(x)A b M) c 



5(i>(x)) 


S(A a M) 




Sfj(x) 


S{ip(y)) _ 


5(A«(x)) 




$v(y) 


5(C a (x)) 


S(A b (y)) 


my) 


SC a {x) 


H( b (y)) _ 


5{Al{x)) 


WW 


5C\y) ' 



0?{x,y) ee (TA a u (x)C b (y)) c = = b£Hr , 
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which, of course, vanish in equilibrium. In these definitions, the symbol T denotes time 
ordering with respect to imaginary time. The identities listed in the right hand sides of 
the previous expressions will be useful later, when we derive the equations of motion for 
2-point functions (see Sec. 3. 3). 

The time ordered propagators can be separated into components which are analytic 



functions of their time arguments [[31], [32]. For example, for the fermion propagator, we 
write 



S(x,y) = 6(t x - T y )S > (x,y) - 9(r y - r x )S < (x,y), 



(2.39) 



with similar definitions for the other 2-point functions D, G, K®, H%, O ab and 0^ b , but 
without the minus sign (remember that the ghosts fields, while anticommuting, are sub- 
mitted to periodic boundary conditions in imaginary time). The two functions S > (x,y) 
and S < (x,y) are analytic functions of x and y for —(3 < lm(x — y ) < and < 
Im (x — yo)</3 respectively. This is the case at equilibrium; we assume that the sources 
satisfy the the conditions necessary to preserve this property off-equilibrium |3T]. We may 
assume, to be specific, that the (analytic) sources are adiabatically switched on at (real) 
time — oo. Then, when (real) time goes to — oo, the various 2-point functions go over to 
their equilibrium values in the absence of sources, and, for arbitrary time, they represent 
correlation functions averaged over the initial equilibrium grand canonical ensemble. 



The induced sources eqs. (|2.23|) - (p.25|) and (|2.33|) - (p.36|) involve products of fields 
with equal time arguments t x = r y . They may be calculated from the time ordered 
propagators (|2.37|) - (|2.38|) taking consistently r y — t x = r\ — ► + . Alternatively, they may 
be expressed in terms of the analytic pieces of these propagators, as defined in ( p.39| ): 



And, 



X 



gYt a K< v (x,x) 



(2.40) 



j^(x) = gTr(^t a S<(x,x)) 
where the trace acts both on spinor and color indices, 



f G M = -9f abc HG<(x,y) 



'cb 



x,y) 



C l a nd {x) = gf abc {& + d£) O bc <(x, y) 



(2.41) 

(2.42) 
(2.43) 
(2.44) 
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Cr{x) = -gr* d»O b ;<(x } y) (2.45) 

where y —* x + stands for r y ~T x — > + . These expressions hold after analytic continuation 
to real time. Then y — > x + should be understood as y§ — Xq = irj and 77 — > + . 

Because the induced sources involve products of fields at the same point, one could 
expect to encounter ultraviolet divergences when calculating them. However, this will 
not be the case here (with the exception of the tadpole term to be discussed in the 
next section). Indeed, as we have seen, in equilibrum all the induced sources vanish, for 
symmetry reasons. Off-equilibrum, no divergence occurs at leading order, as we shall 
verify later. 

For further reference, we list here expressions for analytic components of the free, 
equilibrum propagators. For the fermions, we have 

S<(x,y) = J j^e- ik <^S<(k), S<(k) = $ Po (k)n(k ) = $A(k), (2.46) 

where po(^) is the spectral function 

p (k) = 2ne(k )S(k 2 ), (2.47) 

n(ko) = l/(exp(/%o) + 1) is the fermion occupation factor, and e(k ) = ko/\ko\. Similarly 
the equilibrium free gluon 2-point function Dq (k) is, in Feynman gauge (A = 1), 

D^(k) = -5 ab 9tlu A(k), (2.48) 

where now A(k) = po(k)N(k ), and N(ko) = l/(exp(/5fc ) — 1) denotes the boson occu- 
pation factor. Finally, for the free ghost propagator, we have 

G< ab (k)=5 ab A(k), (2.49) 

because of the periodic boundary conditions imposed on the ghost fields. We note that in 
the description of the plasma using a covariant gauge, the equilibrum state contains ghost 
particles obeying Bose statistics. These are unphysical degrees of freedom, whose role 
is to cancel, when computing physical observables, the contribution of the longitudinal 
gluons, which are also present in the state, according to the propagator ( |2.48j ). 
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3 Equations for the induced sources 



Up to now all the equations that we have written are exact, but incomplete. Indeed, 
these are equations of motion for the mean fields which depend on induced sources given 
themselves by 2-point functions. (A possible contribution from a 3-point function has been 
seen to be negligible in the previous section). These 2-point functions are determined by 
equations which involve generally 3 and 4-point functions, and so on. The approximations 
to be presented in this section will provide us with equations for the 2-point functions 
which involve the mean fields, but no higher point functions, thereby closing the hierarchy. 

The various features of our approximation scheme are summarized in subsection 
3.1. Then, in the following subsections, we go on and derive systematically the equations 
of motion for the various 2-point functions of interest. This is done in several steps. In 
subsection 3.3 we take into account the smallness of the coupling in order to truncate the 
hierarchy. In the next subsections, we implement, for each 2-point function of interest, the 
approximations which follow from the fact that both the wavenumber and the amplitude 
of the mean field oscillations are also small. The technical problems arising from the use 
of a general covariant gauge are discussed in Appendices A and B. 



3.1 Approximation scheme 

The cornerstone of our approximation scheme is the fact that the long wavelength excita- 
tions that we are studying are collective excitations which develop at a particular energy 
scale, gT, well separated, when g <C 1, from the typical energies of plasma particles, ~ T. 
This scale controls both the wavelength and the amplitude of the mean field oscillations. 
Collectivity arises because the wavelength of the oscillations, ~ 1/gT, is much larger than 
the typical interparticle distance, ~ 1/T. 



The way how the scale gT emerges, is perhaps best seen on the equations fl2.18| )- 



( |2.21|) for the mean fields. Technically speaking, one expects collective excitations to 
develop whenever the induced sources on the right hand sides of these equations are 
of the same order of magnitude as the first order differential operators which stand on 
their l.h.s. Assuming that the wavelength/period of the field oscillations is A (i.e., the 
fields carry a typical momentum P ~ gT), the l.h.s. are then of order ~ A _1 0, where 
denotes one of the average fields, that is either (F^) or (ip). The r.h.s., that is, the 
induced sources are indeed of order A~ x when A -1 is of order gT. For instance, eq. (|2 . 2 7| ) 



shows that r\ md ~ S(^). The dominant contribution to S(P) at high temperature is 
~ g 2 T 2 /P. Thus, rf nd ~ (g 2 T 2 /P)(^j), which is ~ P(i/>) ~ A -1 ^) when P ~ gT. 



13 



Similarly, j ind ~ U(A) ~ g 2 T 2 (A) at high temperature, while d(F) ~ d 2 (A) ~ P 2 (A); 
again, j md ~ (9(F) if P ~ gT. Furthermore, the requirement that the plasma stay close 
to equilibrium puts constraints on the strength of the mean fields which also depend on 
A. The existence of such constraints can be anticipated on the basis of the following 
heuristic argument. The plasma will be weakly perturbed if, for instance, the color force 
~ g(F) acting on a typical particle with momentum k ~ T, produces a small change 
in its momentum. Assuming that the force is acting during a time At ~ A, we require 
Ak = g(F)\ <C k. One sees that, for A -1 ~ gT and Ak ^ gT, we have (F) ^ gT 2 . 
More generally, for A -1 ~ gT, there exists a limit in which all terms in the equations 
of motion are of the same order of magnitude. This is achieved when (F) ~ gT 2 and 
~ gT 3 , that is, when these fields are g times a power of T equal to their canonical 
dimension. When the fields take this particular limiting strength, the coupling constant 
disappears from the equations of motion if one expresses all energies and momenta in 
units of gT, and length and time in units of 1/gT. Furthermore, in this limit, all terms in 
the field strength tensor, or in covariant derivatives, are of the same order of magnitude 
(<9 ~ g(A), d(A) ~ g(A) 2 ), and the same is true for all the contributions to the induced 
sources which are non linear in (A). In the rest of this work, we shall work in this limit, 
where the equations keep their full non-abelian structure. However, one should emphasize 
that the equations to be derived are also valid in the case of weaker fields. In the limit of 
weak fields, the equations for (A) and (ip) decouple, and (A) behaves as an abelian gauge 
field, that is, the equations are linear. 

We turn now to the 2-point functions which are needed for the calculation of the 
induced sources and explain the general structure of our approximation scheme. For the 
present discussion, we shall consider as a generic example the off-equilibrium fermion 
2-point function S(x,y). Most arguments extend to other 2-point functions as well. 

The equations of motion for the 2-point functions will be obtained in Sec. 3. 3, by 
differentiation of the mean field equations with respect to the external sources. In doing 
this, we make use of the identities listed in eqs. ( p. 37| ) — ( |2 .38| ) . For S(x,y) one obtains (in 
the case (ip) = ($) = 0) 



S, 1 +g(A) +£) 3 = 1, (3.1) 

where operator notations have been used, and Sq (x, y) = —i(j) x 8(x — y). The off- 
equilibrium self-energy S arises from differentiating the induced source. We argue now 
that such a term can, in fact, be ignored in leading order. To understand why it is so, 
let us recall that we need S(x,y) only to evaluate the induced current jf , according to 
eq. ( |2.41| ). Since the equilibrium propagator S eq (x — y) does not contribute to jf, we need, 
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in fact, only the deviation induced by the gauge fields, 6S = S — S eq . Furthermore, as 
will be clear shortly, we only need 8S(x,y) for small relative distances \x — y\ ^ 1/T, or, 
equivalently, for hard momenta k ~ T. The diagrams contributing to SS involve at least 
one mean field insertion, and the lowest order ones are shown in Figs. 3a and 3b. An ele- 
mentary analysis reveals that, indeed, those in Fig.3b, which represent typical corrections 
coming from S, are negligible compared to that of Fig. 3a when k ~ T. Consequently, we 
shall neglect £ in eq. ( p.l|) and, at the same time, replace S eq by its free value Sq. We 
may then write the final equation as 

(So" 1 +g(4)) S = 1, S = S + 5S. (3.2) 

More generally, in deriving the equations of motion for the 2-point functions in 
Sec. 3. 2, we shall ignore all terms coming from differentiating the induced sources, to 
which we refer as radiative corrections, and at same time replace the equilibrium propa- 
gators by their free expressions. This approximation, equivalent to the familiar one-loop 
approximation to the induced currents (see Figs. 7 and 8 below), closes the hierarchy of 
equations at the level of 2-point functions. The net result is that only mean field insertions 
are to be considered in the off equilibrium 2-point functions. In physical terms, this is 
to say that the dominant interactions in the plasma are those between the hard particles 
and the mean fields. 

It should be stressed that the one-loop approximation alone is not by itself a consis- 
tent approximation for the calculation of the induced sources. It does contain all leading 
terms in the weak coupling expansion, but also terms which are subleading, and which 
should therefore be discarded. These subleading contributions will be eliminated after 
taking carefully into account the other sources of g dependence, namely, the soft charac- 
ter of the average fields. 

In order to implement the approximation that the fields are slowly varying in space 
and time, we need to consider the 2-point functions in real time. As explained in section 
1, these functions, as well as the equations that they satisfy, can be obtained by ana- 
lytic continuation of the imaginary time functions and the corresponding equations. We 
can then proceed to a gradient expansion, for which it is convenient to use the Wigner 
transforms of the various 2-point functions. In the case of the fermion propagator, this is 

S(k, X)=f dh e ik - s S<(X + lx - S -), (3.3) 

where 

T ^ _|_ 

s^ = x^-y^ X^ = — (3.4) 
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Note that in order to avoid the proliferation of symbols, we use similar symbols to denote 
different functions. Thus, for example, we shall often write S < (s,X) in place of S < (X + 
s/2, X — s/2) = S < (x,y). Which function is meant can be recognized from its arguments. 
Wigner functions will be denoted consistently with calligraphic capitals. In equilibrium, 
the non vanishing two point functions are independent of X^ and sharply peaked about 
s M = 0, the range of variation being fixed by the thermal wavelength ~ 1/T. In the 
presence of a slowly varying disturbance, with a wavelength of order 1/gT, we expect the 
s M dependence to be little affected, while the acquired X dependence will have a long 
wavelength variation. This will enable us, for example, to neglect the derivatives with 
respect to X M , of order gT, relative to those with respect to s M , of order T. We shall thus 
transform the equations of motion for the 2-point functions into equations for 5S(k,X), 
valid for k ~ T and dx ~ gT. 

One can use the Wigner function (|3.3| ) to write the induced current ( |2.41|) as 

3?,(X) = gJ-0^Trj»t a 6S(k,X), (3.5) 

where 5S(k, X) = S(k, X) — S eq {k). This expression explains why we need 5S(k, X) only 
for k ~ T: soft momenta k ~ gT represent a small integration region and, consequently, 
give generally a small contribution to such an integral. We are assuming that SS(k, X) 
is not infrared singular, which is the case in leading order. Note that 5S(k,X) is not 
accurately determined by our approximate equation ( |3.2| ) if k ~ gT (then the one loop 
diagrams in Fig.3b contribute at the same order as the tree-level diagram in Fig.3a). 
However, (|3.2| ) provides a SS(k, X) with a smooth behaviour at small k and the integration 



in (3J5) can be extended to all momenta without harm. 

After solving the equations of motion for the 2-point functions, we find that, for 
fields as strong as allowed, the induced pieces of the propagators 5S, 5D and SG, are, 
at most, g times the corresponding equilibrium propagators. It is in this sense that the 
induced effects are perturbatively small, and that the plasma is weakly disturbed. When 
isolating the leading order terms in our equations, we are led to retain terms of lowest 
order in the fermionic fields, that is, terms linear in (if)) and (if)) for K, and 7i, and, owing 
to baryon number conservation, bilinear expressions like (ip)(ip) for 5S and 5T>. However, 
in order to satisfy gauge covariance, we have to keep terms to all orders in (A). In fact, 
the gauge fields enter our equations in two different ways. This may be seen by rewriting 
eq. (|3.1|) as 



+g(A)) 5S = -g(A)S . (3.6) 
The gauge field enters linearly on the right hand side, as a usual perturbation. It also 
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enters non linearly, as part of a covariant derivative, on the left hand side. The existence 
of these covariant derivatives, in the equations for both the mean fields and the 2-point 
functions, will be seen to be the only source of non linearity in (A). This particular way 
in which the mean fields enter our equations is the origin of the fact that there are hard 
thermal loops with any number of gluon lines, but no more than two fermion lines; and 
also of the fact that, in QED, there is no hard thermal loop with more than two photon 
external lines. 



3.2 The induced sources and their densities 



According to eq. (|3.5|) , the induced current jf (X) can be written as an integral over k of 
a density J'^k, X): 

j^(X)=gj ^^(U), J t l(k 1 X)=Try lt f6S(k,X). (3.7) 

Once all 2-point functions are expressed in terms of Wigner functions, this can be gen- 
eralized to all induced sources. In order to facilitate the reading of this paper, we now 
specify the systematics of the notation used for the various densities of induced sources. 
All these densities will be denoted by calligraphic capitals, as the Wigner functions to 
which they are related. 

According to eq. (|2.40|) , the induced source rf nd (X) is given by 



V ind (X)=gJ-0^mX), (3.8) 

where tjL(k,X) = t a 7 M /C^(fc, X). The densities of the induced currents will carry, besides 
Minkowski and color indices, subscripts (f , g, G, Q) and upperscripts (A, tp). As indicated 
after eq. ( |2.32| ), the subscripts refer to which hard particles are contributing to the induced 



sources, f standing for fermions, g for gluons, G for ghosts; Q refers to the tadpole contri- 
bution. The upperscript indicates the nature of the background field which induces the 
current (A for a gauge field, and ip for a fermionic field). We shall see that, in the present 
approximation, one can separate the contributions of the fermionic and gauge fields, (up 
to trivial color "rotations" imposed by the gauge transformation properties of the various 
densities). Thus, we shall have a collection of induced currents {jf, j£, j£, j£, jf , j^}, 
together with the corresponding densities {^T f A , j7 g A , Jq , J^, jf ', J%}- In some cases, 
it will be convenient to operate various grouping of these contributions. The absence of 
the upperscript will imply a grouping of all contributions corresponding to the various 
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upperscripts: 

Ji — J A + Jf , Jg = Jg+Jg, Jg = Jg i Jn = Jn- (3.9) 

Similarly, the absence of a subscript implies a grouping of all contributions corresponding 
to the various possible subscripts: 

J A = Jf A + J* + j£ + j£, J^ = 3t + Ji- (3.10) 

We shall also use the subscript b to group together all the contributions coming from 
bosonic particles: 

Jt = Jgi = J? + Jg+ Jn, (3.H) 

and similarly for the corresponding induced currents. 

We shall also be led to consider the baryonic vector and axial currents. They are 
defined by 

b^x) = Tr^5S<(x,x) : b%x) = Tr l5lll 6S < (x,x) (3.12) 

and their densities will be denoted by B^(k,X) and B^(k, X), respectively: 

B,(k,X) = Tr lfM 5S(k,X), B%k,X) = Tr l5lfl 5S(k,X). (3.13) 

Finally, as we shall see, a fermionic field may also induce oscillations in the axial color 
currents jf and jj, with 

j?*(x) = 9 Tr ( l5l ,t a 5S<(x,x)) , J { 5 ;(k,X) = Tr { l5ltl t a SS{k, X)) . (3.14) 

and similarly for j^. 

We end this subsection by calculating the induced current jo which arises from 
the four gluon interaction. According to eq. (|2.36|) , jn is proportional to the gluon loop 
{A a A b v ) c = D°^,(x,x + ), where D^ u (x,y) denotes the (exact) gluon propagator. Since this 
is multiplied by g 2 , in order to get the leading order contribution to jn we can replace 



it whith the free, equilibrum propagator, given by eq. ( 2.48 ). Indeed, according to the 
previous discussion, the difference D — D is expected to be, at most, of order gD . We 
obtain then 

3W = g 2 tf v %{A^x)) Dfc d {x,x+) = -fi 2 (A;(x)), (3.15) 
where we have noted 
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and eq. Q2.48| ) for Dq°£ in Feynman gauge has been used. (Note that the ultraviolet 
divergence associated with the T = part of D (x,x + ) has been removed by the usual 
zero-temperature substraction.) The momentum space density for ja is 

J&Jk,X)= -gN(d-l)A(k)(A a (X)). (3.17) 



3.3 Equations for the 2-point functions 

In this section, the imaginary time equations of motion for the 2-point Green's functions 
are obtained by differentiating the mean field equations fl2.18| )- (|2.21 ) with respect to 



the external sources. As explained in the previous subsection, we ignore all radiative 
corrections, that is, all terms coming from differentiating the induced sources. 

We define 



g^di-{l-\)(f^\ 5 ab 

+gf abc [r, upX (A?(x))d x x - r, pXv (dZ(A&x)))] - g 2 t^%(A^x))(A' d (x)). (3.18) 



This operator appears whenever we differentiate eq. (|2.18|) with respect to an arbitrary 
source J 

wr+^A(ffl)} = (3 - 19) 

and it may be interpreted as the inverse propagator of a gluon in a background gauge 
field (see eq. ( ET27| ) below). 

In writing the equations of motion for the 2-point functions, we shall ignore, when- 
ever they appear, the following abnormal propagators (Tijj(x)ip(y)) c , (Tip(x)((y)) c , 
(Tip(x)((y)) c , together with their Hermitian conjugates. These functions are vanishing in 
equilibrum, while their off-equilibrum expressions involve at least one pair of weak fields; 
because of this, they are found to be negligible relative to the leading terms appearing 
next to them in the equations. In order to make this statement clear, it is useful to analyze 
a specific example. Consider the equation for K®(x,y) which follows from differentiating 
eq. (|2.18|) with respect to fj 

C%{y)K v h {x,y) = gS(x,y) % t a {^(y)) + g^{y)) ltl t a (T^{y)^{x)) c 

-gf abc (d»(( b (y))) W(x)C c (y)> c - gf abc (d^(T^(x)C b (y)) c ) (C(y)). (3.20) 

The first term in the r.h.s. of this equation involves the off-equilibrum quark propagator 
S(x,y). We write it as S = S + SS, where So is the free equilibrum propagator, while 



19 



<55' is the correction due to the average fields. We assume that 5S ^ gS , and this will be 
verified explicitly in Sec. 3. 4. In leading order, we can therefore neglect the contribution 
of SS relative to that of So- As mentioned, the other 2-point functions in the r.h.s. of 
eq. ( p.20| ) are also expected to be small with respect to So. We show how to verify this 
explicitly for the first of them, which is proportional to (Tip(x)ip(y)) c . By differentiating 
eq. (|2.19|) with respect to f/, we obtain the equation 

% (T^{x)^{y)) c = igK a u (y,x)Yt a (Hx)). (3.21) 

As we shall see in Sec. 3. 4 (after eq. ( |3.47| )), this equation implies that (Ttp(x)ip(y)) c ~ 
(So/T 3 )(ip) 2 . Therefore, the second term in the r.h.s. of eq. ( |3.20| ) is g{^(y)) r ) l it a 
(Tip(y)ij;(x)) c ~ gS ((ij)(ij)/T 3 )(tJj). Relative to the first term gS (x - y)7^ a (?/%)), 
this is of order /T 3 , that is, at most of order g. Similar considerations apply 

to the last two terms in eq. ( |3.20| ), which can be also neglected. Our final equation for 
K%(x,y) is therefore 

C%{y)KZ{x,y) = gS (x - y)^t a (^(y)). (3.22) 

The previous arguments may be extended to all the other equations of interest. We 
shall not repeat them, but simply list the equations obtained after the corresponding 
simplifications have been done. 

By differentiating again eq. Q2.19D , but this time with respect to j b , we derive a 
second equation for K b (x, y) 

%Kl(x,y) = -igt a j^(x))DZ b {x-y). (3.23) 

A diagrammatic interpretation of eqs. ( 3.22| ) and (|3.23 ) is given in Fig. 4. There are, of 
course, similar equations for but we do not write them explicitly, because they are 
simply related (by Hermitian conjugation) to eqs. (|3.22|) and ( |3.23j ). 

The equation for S(x, y) is obtained by differentiating eq. (|2.19| ) with respect to r\ 
% S(x, y) = i5(T x - r,)5 3 (x - y) - igj v t a Oj(x))H^x, y). (3.24) 

The diagrammatic contributions to 5* are displayed in Fig.5. 

By differentiating the equations Q2.20| ) and ( [2.2 1|) for the ghost fields with respect 
C b and C a , we generate two equations for G ab (x,y): 

dl G ab (x, y) - gr d (A^x)) d»G db (x, y) = 

S ab 5(r x - r,)5 3 (x - y) - g r d d» \(( c (x)) Of(x,y)} , (3.25) 
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and 



% G ab (x, y) - gf bcd {Al(y)) d»G ad (x, y) = 

S ab 5(r x - r,)5 3 (x - y) - gf bcd (><C%)>) O ad (x, y). (3.26) 



Finally, by differentiating eq. (|2.f8| ) with respect to we get the following equation 



for H*{x, V ): 



fJ,U 

ab 



C%(x)D2(x,y) = g,u5 ab 5(r x -r y )5\-K-y) 

+ g{^(x))^t a K b {x,y)+H b (y,x) 1 ^ a (i>(x))} 

- gf acd {(d,(C(x))) O db (x,y) + (d;O b v c (y,x))(C d (x))},(3.27) 

Diagrammatic contributions to D are shown in Fig. 6. 

We have not written here the corresponding equations for the ghost-gluon abnormal 
propagators, O ab and Of. We do this in Appendix A, where we prove that, in leading 
order, the equation of motion for the ghost mean fields coincide with the tree level equa- 
tions; there is no collective excitation with ghost quantum numbers. Accordingly, we shall 
not be concerned anymore in the main text with the propagation of the ghost fields. We 
set the ghost sources C a and C a to zero. Then (() and (Q vanish, and all terms involving 
O ab and Of disappear from eqs. (gg§-(gg7[) . 

In what follows, we simplify the writing and denote the average gauge potential or 
the average fermionic field respectively by A(x) and ip(x), omitting the brackets. We also 
denote simply by D or F^ v the covariant derivative or the field strength tensor associated 
to the average gauge field, in place of T> or T^. Calligraphic capitals will therefore be 
reserved for Wigner functions and for the densities of the induced sources. Furthermore, 
we shall restrict our discussion in the main text to the Feynman gauge (A = 1); we prove 
in Appendix B that our final equations are independent of the gauge fixing parameter. 

The calculation of the induced sources from the 2-point functions given by the 
previous equations is equivalent to the calculation of the one-loop Feynman diagrams 
displayed in Figs. 7 and 8. The approximations to be implemented in the rest of this 
section will provide us with consistent equations from which the leading contributions to 
these diagrams will be directly obtained. 



3.4 Equations for /C and TL 

The real-time equations which determine the abnormal quark-gluon propagator K£ are 
obtained from the imaginary-time equations ( |3.22| ) and ( |3.23| ). We consider r y > t x and 
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obtain, after analytic continuation, 



CZ(y) K» b <(x,y) = -gS<(x - yh,t a ^(y), 



(3.2* 



and 



p x K^<(x, y) = -igt a ^(x) D^(x - y). 



(3.29) 



By acting on eq. ( |3.29|) with the differential operator p x from the left, we obtain 
p x p x K: < (x,y)=igp x [t a Y^(x)A(x-y)}. (3.30) 
To simplify the l.h.s. of this equation, we use the identity 



p x P, 

where o» v = {i/2)[^,Y} and 



-a 



Ffiuix), 



(3.31) 



D 2 = d 2 + 2ig A(x) ■ d x + ig(d ■ A(x)) - g 2 A 2 (x). 
Then we replace x and y by the coordinates s and X (see eq. 



(3.32) 
and rewrite the 

derivatives as d x = d s + \dx, d y = —d s + \d x and d 2 x — d 2 = 2d s ■ d X - We recall that 
a derivative dx is to be regarded as being of order gT, and a derivative d s as being of 
order T. Furthermore, the mean field is assumed to be of order T, that is gA is to be 
considered of the same order as the derivative dx, which insures the consistency of the 
soft covariant derivative, Dx = dx + igA. As for gF^, it is of order g 2 T 2 . Finally, we 
expand A(x) = A(X + s/2) ~ A(X) + (s/2) ■ dxA(X), where the second term is g times 
the first one, and likewise for the other fields. Finally, eq. (|3.30| ) reduces to 



~d 2 + d s -d x + 2igA(X) ■ d s ] KZ<(s,X) = ^ S A) 7 ^>(X). 

Similar manipulations performed on eq. ( |3.28 ) lead to 

5 a V - 9 S ■ d x ) - gf abc T^ pX AP(X)d*} K»<(s,X) 

= -gS<(s) llM t a if;(X). 



(3.33) 



(3.34) 



We remark that the contribution from the four gluon interaction (i.e., the term propor- 
tional to f ) disappears at this order. 

In writing eqs. ( |3.33|) and (|3.34j) , we have neglected terms of order g 2 T 2 inside the 
brackets on the left hand sides. As for the remaining terms, the first one, d 2 , is of order 
T 2 , while all the others are of order gT 2 . Therefore, in order to have consistent equations 
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order by order in g, we should also expand K< = K<® + + with K« ~ gK^\ 
etc... Of course, because of the approximations which have been done, eqs. ( |3.33|) and 
( |3. 34 ) can only be used to determine . The dominant terms in the two equations lead 
to the consistency condition 



d 2 K^°\s,X) = 0. 



(3.35) 



Next, we note that the "big" derivatives d 2 ~ T 2 give rise to contributions of the form 
d 2 K^ which are of the same order as, for example, (d s ■ dx)K^\ Such terms are incon- 
sistent and, besides, they cannot be easily evaluated. In order to do so, one would need 
equations for K^ 1 ' containing the unknown function K^ 2 \ and so forth. However, such 
"spurious" contributions cancel in the difference of the two equations ( |3.33| ) and (|3.34| ). 
By doing this difference, we end up with the following equation for K^°\ which we denote 
from now on simply by K < : 



5 a V (dx + igA(X)) ■ d s + g -r bc T, vpX A" c {X)dl 



s,X) 



^[S<(s)+0 s A(s)h,mX). 



(3.36) 



When Fourier transformed with respect to s, eq. (|3.36j ) translates into an equation for the 
Wigner function /C^(fc,X): 



k-(d x + igA(X)) K.%(k,X) + p ab T, upX A^X)k x mk,X) 
= i 9 -p G {k)[N{h)+n{h)}$ lp t a <iP{X). 



(3.37) 



From the right hand side of eq. ( |3.37| ), one sees that JC® is proportional to fc5(k ). Con- 
sequently 



k 2 !C*(k,X) = 0, 



(3.38) 



The first identity above is the consistency condition ( 3.35 ). It simply expresses the mass- 
shell condition for Kft. The second identity is the Dirac equation for the quark-gluon 
propagator. It can be also verified that the right hand side of eq. (|3.37| ) is transverse with 
respect to k, so that 



k»lC a Jk,X) = 0. 



(3.39) 



By using these properties, namely, Eqs. ( [3. 38| ) and ( |3.39|) , one easily transforms eq. (|3.37f) 
into an equation for Jfi a (k,X) = ■y tl IC p (k, X) 

k ■ \{d x + igA{X)) 5 ab + igA ab {X)] l£ b (k, X) 



ig^p (k)[N(k )+n{k m a ij;(X), 



(3.40) 
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where the factor (d—2) results from the identity 7^7^ = — (d— 2)}k. Eq. ( |3.40| ) is covariant 
with respect to local gauge transformations of the mean fields ip and A. This is seen by 
noting that the differential operator from the left hand side is a covariant derivative acting 
on vectors in the direct product of the fundamental and the adjoint representations. It 
follows that ¥fi a transforms as t a tp(x) under a gauge transformation. Explicitly, when 
ip(x) — > S(x)ijj(x), with S(x) = exp(i9 a (x)t a ), then 

ty°(k,X) - ~S ab {X)S l3 {X%){k 1 X) 1 (3.41) 

where S(x) = exp(i6 a (x)T a ) is the gauge transformation in the adjoint representation. 

The equation satisfied by tfi(k,X) = t a fjt a (k,X) results from eq. ( |3.40|) after multi- 
plication by t a 

{k ■ D x ) %{k, X) = -ig^zlc fPo (k)[N(k ) + n{k )\mX), (3-42) 

where Cf is the quark Casimir. The covariant derivative which remains in the left hand 
side of eq. ( p.42j ) is that of the fundamental representation, Dx, showing that If, transforms 
as ip. In contrast with what happens in eq. ( |3.4(J| ) for }f a , the gauge field insertions on the 
hard gluon lines have disappeared in eq. ( |3.42| ), which coincides with the corresponding 
equation for an abelian plasma [[If] . 



We end this section by writing the equation satisfied by the Wigner transform of 
the quark-gluon propagator H^ K . This follows from eq. ( |3.42j ), by noting that Tflik, X) = 
^(fc,X) 7 °, 

T/C(k, X) (k ■ D x ) = ig ^-lc fPo (k)[N(k ) + nfaWXW. (3.43) 
Here = — igA^, where the derivative operator = ( d/dx^j acts on its left. 



3.5 Equations for 5Sa and 5S, 



The equations for the off-equilibrium fermion propagator are obtained by analytic con- 



tinuation from eq. (|3.24j) and its conjugate: 

p x S<(x,y) = ^Wffrfxj), 



(3.44) 



S<(x,y)pl = -igK<(x,y)i>(y)t a Y, (3.45) 

with Iffi = ffi — igA and (jfi = (^d/dx^ 7 M . We transform them by first acting with 
the covariant derivatives Jj) x on eq. ( |3.44|) from the left, and with Jfr y on eq. ( |3.45| ) from 
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the right. We subtract the resulting equations and separate the leading order terms by 
manipulations similar to those encountered in the previous subsection. A new feature 
stems from the fact that in equilibrium S < (x,y) does not vanish, but is equal to Sq, 



eq. (|2.46|) . As explained in Sec. 3.1, we write S < (x, y) = Sq (x — y) + SS(x, y) and assume 



SS to be of order g Sq. It is then straightforward to derive the equation satisfied in leading 
order by the Wigner function 5S(k,X), 



[k-D x , 5S(k,X)) = -i|FH<V, S<(k)]-g(d x A^) k,d k u S<(k) 



i 9 - {K:(k,X)^(X)t a Y^ - Wt a i>{X)H a v {k,X)} . 



(3.46) 



We call 5S (k, X) the solution of this equation when ip = ip = 0, and SS^(k,X) 
the correction which includes the effects of the fermionic fields. We have 5S(k,X) = 
5S A {k,X) + 5S^{k,X), and 

k-D x , ss*{k,x)] = -t^{}ci(k,x)^(x)t a Y^-^ u ty(x)n:(k,x)}. (3.47) 



By using this equation, together with eqs. ( |3.42|) and ( |3.43|) , we can verify that 5S^(k, X) ~ 
(So/T 3 )ip ip. Thus SS^(k,X) < gS when xjjijj < gT 3 , as anticipated in Sec. 3.1. Remark 
also the similarity between eqs. ( |3.44|) — ( 3.45|) for S < (x,y) and eq. ( |3.21| ) for the abnormal 
propagator (Tip(x)ip(y)) c . By following the same steps as above, one can derive for the 
Wigner transform of (Ti[)(x)i/j(y)) c an equation analogous to eq. (|3.47|) for SS^(k, X). 
From it we derive the estimate (ipip) c ~ {So/T 3 )ip 2 which has been used in Sec. 3. 3 in 
order to discard the contribution of such abnormal propagators. 

As a color matrix acting in the fundamental representation, SS^ can be decomposed 
into a color singlet and a color octet: SS^ = SS^ + SSft a . Besides, SS"^ is a 4 x 4 matrix 
in Dirac indices. We can verify on eq. ( 3.47 ) that trSS^ = ti'y^SS^ = tia^dS^ = 0, and 
these equalities hold for both the color singlet and the color octet (the trace acts on Dirac 
indices only). The quantities tvy ^SS^ and tvy^'j ^SS^ are nonvanishing. Thus a fermionic 
field induces oscillations in the vector and axial baryonic currents and 6^, eq. (|3.12|) , 
as well as in the vector and axial color currents, jf and The corresponding densities 
are given in eqs. ( |3.13| ) and (|3.14| ), with 5S replaced by 5S^(k, X). 

Let us consider now the equation satisfied by 8S A {k,X). By using eq. (|2.46|) , we 
can write it as follows: 



k ■ D x , SS A (k, X) =-g{[k- D x , Aj + ftk, {d x A») d u } A(k). 



(3.48) 



This equation shows that 5S A ~ g(A/T)So, so that SS A ~ gSo when A ~ T. It can also 
be verified that 5S A oc 7 M t a (Tr^t a 5S A ). Thus, a gauge field induces only fluctuations in 
the vector color current. 
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Eq. ( P-48| ) is not covariant under local gauge transformations. This reflects the lack 
of covariance of our definition for the fermionic Wigner function 5S(k,X). We could, 
instead, define a gauge covariant Wigner transform, as for example in ref. j2"2f . But it is 
easier to recover gauge covariance directly at the level of eq. ( |3.48|) . This can be done in 
a standard way, by replacing the canonical momentum k, associated with the gradient 
operator, by the kinetic momentum p = k — gA, related to the particle velocity. By 
applying this operation to the quark Wigner transform S(k, X), we obtain a new function 
of p and X, <S(p,X) = S(k = p + gA, X), which transforms covariantly. This can be 
verified directly by writing (for ip = ip = 0) S(p,X) = Sq(p) + 5S A (p,X). Then, in 
leading order, 

5S A (p,X) = 6S A (p,X)+g(A(X) ■ d p )S<(p), (3.49) 
and 8S A {p, X) satisfies the gauge covariant equation 

p-D x , 5S A (p,X)] =gp"p-F(X)-d p A(p), (3.50) 



deduced from eq. ( [3.48 ) . This is the non-abelian generalization of the well-known lin- 



earized Vlasov equation. Note that although gA ~ gT <C k ~ T, the change from k 
to p in S(k,X) is non trivial, as the second term in eq. ( |3.49| ) is of the same order of 
magnitude as the first one. However, this term is a total derivative with respect to p, so 
that it does not contribute to the induced current jf (see eq. (|3.5| )). 

The change from canonical to kinetic momentum in SS^(k,X) is irrelevant at this 
order. Indeed, if we replace k by p = k — gA in eq. ( |3.47|) , we only generate terms of 
higher order. And, in fact, eq. (|3.47|) is gauge covariant as it stands. 

In summary, the previous discussion shows that all effects induced by the mean 
fields on the hard quarks can be described by the gauge covariant Wigner function 

S(k,X) = S<(k) + 5St(k } X) + (SS*(k,X) +5S A (k,X)) t a , (3.51) 

and k is to be understood here and from now on as the kinetic momentum. In particular, 
the induced current jf has the gauge covariant density 

Jt%(k,X) = Tr lp t a (NrfS^X) + 5S^(k,X)) = J Aa (k,X) + J?;(k,X). (3.52) 

The factor Nf reflects the fact that all quark flavors contribute identically to the current 
induced by the average gauge fields. Similarly, eq. ( 3.47 ) contained an implicit sum over 
the different flavors of the quark mean fields. 
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3.6 Equations for 8Q 



From eqs. ( |3.25| ) and fl3.26| ), one obtains the following real-time equations of motion for 
the ghost propagator 



and 



d 2 x G<(x, y) = gr d A c (x) d£G%(x, y), 



d 2 v G< b (x, y) = gf bcd A c (y) d»G< d (x, y). 



(3.53) 



(3.54) 



After separating the induced part, by setting G^ b (s, X) = GQ ab (s) + 5G ab (s, X), with 
Go a b given by eq. (|2.49| ), we transform these equations in a way which should by now be 
familiar. We obtain 



k ■ 



d x 5g ab (k,X) + 9 - (f acd 5g cb + f bcd 6Q ac ) A d (X) 



2 

9 rabc 



f abc (d u A^(X))k^ k A), 



(3.55) 



where SQ(k, X) is the Wigner transform of 5G(s,X). This equation shows that, as ex- 
pected, SO ~ g(A/T)A. 

For later reference, we note that 8Q ab determined by eq. (|3.55|) is a color matrix 
of the adjoint representation, i.e., it can be written as 5Q ab = (T c ) ab 5Q c , where the new 
functions 5Q a (k,X) satisfy 

(k ■ d x )5g a (k,X) - 9 -f abc (k ■ A b (X))5G c (k,X) = -| (cU*(X)) k^A). (3.56) 

Then the ghost contribution to the induced current density may be written as (see 
eq. (HD) 



j£Jk,x) = -ik,r bc 5g bc (k } x) = -Nk,sg a {k } x). 



(3.57) 



3.7 Equations for 5T>a and ST>^ 

After analytic continuation to real time, eq. (|3.27|) gives 

C™(x)D%<(x,y) = g {$(x)^t a K b u <(x,y) + H b > (y,x)^(x)} . (3.58) 

Another equation for D K is obtained by interchanging in this equation the space-time 
arguments, as well as the color and Minkowski indices, 

C b ;Ay)DZ: < (x 1 y) = g{^(yha b K; > (y,x) + HZ<(x,y)^ (3.59) 
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The symmetry property D°%, < (x,y) = D b Z > (y,x) has been used. We set D°% / < (x,y) = 
—g flv 8 ah A(s) + 5D"^,(s,X), with 5D expected to be of order gA, and we separate the 
leading order terms in eqs. ( |3.58| )-( |3.59| ) in order to get equations for the Wigner function 
5V p b u (k, X). In order to do this, we need all the analytic pieces of the abnormal quark- 
gluon propagators, i.e., not only K < and H < , for which we have obtained equations in 
Sec. 3. 3, but also K > and H > . In fact, by repeating the calculations of Sec. 3. 3, one can 
show that, in the present approximation, K > and K < satisfy identical equations, and 
both vanish in equilibrium, so that 

K$>(k,X) « K?<{k,X) = K*(k,X), (3.60) 

and similarly for the 7i functions. 

In order to separate the effects of gauge fields from those of fermionic fields, we 
write 8T> = 5T> A + 5T>^, where 5D A (k, X) is the modification induced by the gauge field 
when if) = = 0. At leading order in g, eqs. (|3.58|) and ( [3.59D imply the consistency 
conditions 

k 2 6V Aa u b (k, X) = igf abc A{k)Y, upX A" a {X)k\ kW; b = 0. (3.61) 

These will be used later in simplifying some formulae. At next to leading order, we take 
the difference of eqs. Q3.58|) and ( |3.59|) and obtain, after a lengthy, but straightforward 
calculation, 

k ■ d x bv A ; b - 1 (r d v, aXp 5v A d : b + f bcd v uaXp 5v A r) & k 

= -i 9 - j ahc r M , P A {8T X A§ k\d k a A), (3.62) 

and 

k ■ d x SV+* - I (f acd T, aXp SD+£ + f bcd T uaXp 61>+Z) k x A p c 

= % g -i, { lp x a K, b v - 7 ,t b /c;} + 1 1 [n b vlfi t a - n; lu t b } (3.63) 

Note that the four gluons vertex does not contribute to this order and in this particular 
gauge, so that there are no terms quadratic in A in the right hand side of eq. ( [3.62] ). 

The contribution of hard gluons to the induced current can be calculated from 
eq. ( |2.43| ), which leads to the density 

J^(k, X) = i f abc T ppXv k x 5V p b :(k, X). (3.64) 

The separation 5T> = 8T> + 5T>^ leads to a corresponding separation in the induced 
current density, J g = J A + . We are interested now in transforming the previous 
equations for 5T> A and 5T>^ into equations for these densities. 
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Let us consider first the density J"^ . This is obtained from SD^, given by eq. ( |3.63| ) 
above. Since the r.h.s. of this equation is antisymmetric with respect to simultaneous 
permutations of the color and Minkovski indices, the same property holds for 5T>^ 



SV^\k,X)= -SV% a (k,X). 



(3.65) 



Furthermore, the right hand side of eq. (|3.63|) is transverse with respect to k. Indeed, by 
taking its scalar product with k^, we find in the r.h.s. combinations like k ■ K, a or ftlC® 
which are vanishing in the present approximation, as shown by eqs. ( |3.38| ) and (|3.39| ). 
This implies 



k»5Vt ab (k,X) = 5Vt ab (k,X)k» = 0. 



(3.66) 



Using this property, together with the definition Q3.64T ) for the current density J^£, as well 
as the explicit form of the T symbol, eq. ( |2.13|) , one obtains 



j^;(k, x) = -i k, r bc (svty (k, x) 



(3.67) 



By using the transversality property Q3.66j ), we readily derive from eq. (|3.63| ) the 
following equation for (SV^ : 

k ■ d x (5Vt b )l - g(k ■ A c ) {f acd (5V*)l + f bcd (8Vt d )l} 
= i 9 -$ {t% b - t h 1f, a } + % | [lfi h t a - lfi a t b ) i). 

From eq. (|3.40|) , one can show that 

t% b - t b }f a = i f ab % c = -T c p. 



(3.68) 



(3.69) 



It then follows from eq. ( |3.68|) that (ST>^ is a color matrix in the adjoint representation, 
that is, we can write 

(5Vt b )l = ~(T c U5Vt = -8V*,, (3.70) 
and thus defined satisfies the gauge covariant equation 

~k-D x , 6V*(k,X)] = i^T a {i;(X)ljC a (k,X)-lfC a (k,X)4>(X)}, (3.71) 

where [£)„ O] = d^O + ig[A m 0} for any 0(X) = O a (X)T a . 

It is now straightforward to write the equation giving the current density . From 
eqs. ( |3.67| ) and ( |3.70| ) we have 

J^ a {k,X) = k^TrT a SV^(k,X) = Nk^V^k.X). (3.72) 
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Then, from eq. ( |3.71| ) we deduce the following gauge covariant equation for J~£ = t a J^'- 



k-D x ,jfJk,X) 



i 9 -Nk^ [${X) IjC a (k, X) - Tfi a (k, X) i/>(X)} • (3.73) 



There is, of course, a similar equation for the density J7g 5 of the axial color current induced 
by fermionic fields on the hard gluons. 

Let us consider now the color current induced by the gauge fields, i.e. j^. Its density, 
J"/, is related to 5V A by eq. From eq. we see that SV^ b is symmetric in 

the Minkovski indices and belongs to the adjoint representation, that is 



(3.74) 



with the new functions SVt v satisfying 



k ■ d x 5V* a - 9 -f abc A p b k x 



ucrXp 



fiupX 



In terms of these functions, the density j7" g A is 

j£(k,x) = N{svfi(k,x)r - k,5vt:(k,x)} 



(3.75) 



(3.76) 



and it depends only on the particular combinations SV^k u and ST>^. The equations 
satisfied by these quantities are easily obtained from eq. (|3.75|) . They read 



and 



(k ■ d x )SVt v - g r bc ik ■ A b )5V c / + g f abc At [bV%k 

= g(d — 1)^ (d x A%) (d*A), 



(k ■ d x ) (8VIX) - 9 f abc (k ■ A b ) (SV^k 1 



(3.77) 



+ y -f abc {k,(A b ■ 6V C ■ k) + Al(k ■ 6V C ■ k) + k?AZSV° vlt + (k ■ A b )k»6V^} 

= kk,K - g,uk 2 ) {d x AD (d k p A). 



(3.78) 



We have dropped the superscript A on 5D in order to simplify the writing, and this will 
be done in the following equations as long as no confusion is possible. Furthermore, we 
have used compact notations for tensor products, such as A b ■ 5V C ■ k = A p h bT> c pv k v . 

We have to solve eq. ( |3.78[) before eq. ( j3.77|) . Indeed, its solution, dV^k", enters 
the l.h.s. of eq. (|3.77|) for 8T>® U . In order to do so, we first replace the term proportional 
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to k 2 5V c ^ in the l.h.s. of eq. ( gTg ) by the expression deduced from eq. ( |3.61| ). Second, 
we note that the r.h.s. of eq. Q3.78 ) is transverse with respect to k, so that 



k t *6T% v {k,X)k u = 0. (3.79) 

This allows us to eliminate the corresponding term in the l.h.s. of eq. ( |3.78j ). We are then 
left with the equation 

(* • d x ) (6V; U r) - \g f ah \k ■ A b ) (5V% k») + 9 -f abc k,(A b ■ 6V C ■ k) 

= -AKK - g, v k 2 ) (d p x A v a ) (SjA) - -a 2 f abc {k ■ A b )A^A. (3.80) 



The structure of this equation suggests the decomposition 8V^ v k u = gA°A + where 
8V a are new functions which satisfy 



(k ■ d x )5V; - \gf ab \k ■ A b )5V^ + 9 -f abc k,{A b • 6V C ) = ^M^xO^A). (3.81) 



By comparing this equation with eq. ( |3.56|) for the ghosts propagator, we see that 8T> a 



kfj,SQ a . Therefore eq. (|3.80|) has the following solution 



6V; v k"(k,X) = gA a u (X)A{k)-k,6g a (k,X), (3.82) 



which, not surprisingly, relates ghost and longitudinal gluon degrees of freedom. We use 
this relation in order to rewrite the induced current density Q3.76Q as 



j£(k, X) = gNA{k)A«(X) - Nk^ (6V^(k, X) + SQ a {k, X)) , (3.83) 



where the superscript A has been reintroduced. The relation (|3.82j) can also be used to 
transform eq. ( |3.77| ) for 8V A " into 



k ■ D x , 5Vr(k,X)\ = g(d - l)Av(M)(^A) + gf abc (k ■ A b )5G c (k,X). (3.84) 

We now consider the density J A of the total color current carried by gluonic degrees 
of freedom, j A = j A + jc + jn- By adding eqs. ( |3.83|) , ( |3.57|) and ( [3.17] ), we obtain 



J h A ;(k,X) = -gN(d - 2)A(k)A;(X) - Nk^ {5V A »(k,X) + 25$ a (k,X)} . (3.85) 

The first piece in the right hand side is the tadpole contribution to the induced current 
density, including only the N(d — 2) transverse gluons. The second piece contains the 
contributions involving the 3-gluon vertex. In this term, the role of the ghosts is clearly 
seen: the term 25Q a compensates for the two spurious degrees of freedom contained in 
5V Av . 

av 
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The equation satisfied by j7~ b A follows from eqs. ( |3.56| ) and fl3.84| ) 

k ■ D x , J h \(k, X)} a = -gN(d -2){[k- D X , A,] a + k,k v (d p x A u a ) d k p } A(k). (3.86) 

This equation is not gauge covariant, in spite of the covariant derivative in the l.h.s. 
and the presence in the r.h.s. of the explicit factor N(d — 2) indicating that only the 
transverse gluons contribute. This lack of covariance has its origin in the noncovariant 
character of the Wigner functions that we have been using. Recall, however, that the 
physically relevant quantity is the induced current j£, which is obtained from its density 
after a k integration. This means that j£ is determined only up to a total derivative with 
respect to k. We use this arbitrariness in order to construct a gauge covariant density. 
We define, in analogy with eq. (|3.49|) , 

J b *(k, X) = J h A ^(k, X) + gN(d -2){A- 8 k ){k^{k)). (3.87) 

This new density gives the same current as J^(k,X) after the fc-integration, and it 
satisfies the gauge covariant equation, derived from eq. ( |3.86| ), 

~k-D x , J^(k,X)] a = gN(d-2)k tl k»F^d»A(k). (3.88) 

We recognize here a situation similar to that encountered in Sec. 3. 5, in relation with 5Sa- 
By analogy, we can interpret eq. ( |3.87|) as expressing the change between canonical and 
kinetic momenta for the N(d — 2) transverse gluons. 

This may be seen more explicitly by expressing the current density = + j£ 
in terms of a gauge covariant Wigner function, denoted by T)(k, X). This is a (iV 2 — 1) x 
(iV 2 — 1) color matrix whose diagonal elements may be interpreted as the distribution 
functions of transverse gluons of each color. In equilibrium, it reduces to 

V^{k) = 5 ab (d - 2) A(k), (3.89) 
since there are (d — 2) polarization states for gluons of each color. Off-equilibrium, 

V(k,X) = T) (k)+5T)(k,X), (3.90) 
where 5V(k, X) is related to the induced current jb = ih + it by 

Jb yX) = gJ-^k,trT a 5V(k,X), (3.91) 



(2nY 

with the trace acting in the color space. Since ST>(k, X) belongs to the adjoint represen- 
tation, it is completely determined by eq. ( |3.91| ). In line with the previous analysis, we 
write 

SV(k,X) = 5V^{k,X)T a + 5Vt{k,X)T a . (3.92) 
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where SV^ satisfies the gauge covariant equation ( |3.71| ), while 5V A obeys an equation 
analogous to that for J A , ( ggp , 



k ■ D x , 5V A (k } X)\ =g(d-2) k p F pu d v A(k), (3.93) 
with F vp = F? p T a . 

In Feynman gauge, the relation between ST> A (k, X) and the gluon and the ghost 
propagators is given by (see eq. ( |3.85|) ) 

5V A (k,X) = -{6Vi"(k,X) + 2 5g(k,X)} + g(d-2)(A(X)-d k )A(k). (3.94) 

If we also notice that, in the same gauge, T>a(k) = —D^ik) — 2Gq (k), we see that 
X>(fc, X) may be viewed as resulting from the replacement of the canonical momentum by 
the kinetic one in 

V(k, X) = - {V»(k, X) + 2Q(k, X)} . (3.95) 

That is, T>(k,X) « T>(k + gA,X) in leading order, and T>(k,X) may be interpreted as 
the distribution function for the transverse gluons. Of course, the separation of T>(k, X) 
in terms of T> v v and Q depends upon the choice of the gauge. However, equations ( |3.73j ) 
and ( |3.88| ) for the induced currents, as well as the corresponding equations ( |3.71| ) and 



( |3.93p for the covariant Wigner functions, are independent of gauge fixing, as shown in 
Appendix B. 



4 Kinetic equations for hard particles 

In the previous section, we have derived equations for the Wigner transforms of the off- 
equilibrum 2-point functions which are needed in the calculation of the induced sources for 
the average fields. These equations describe the dynamics of the plasma particles in the 
presence of the average fields which may be considered at this stage as given background 
fields. We analyze now the structure of these equations and we write their solutions for 
retarded conditions, that is for average fields which vanish adiabatically as Xq — > — oo. 
As we shall see, the equations for the densities of induced sources take the form of kinetic 
equations, that is, of equations for generalized phase space distribution functions for the 
plasma particles. 
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4.1 Equations for the densities of the induced sources 



The equations for the densities of the induced sources have been already written in the 
previous section, or can be derived in a straightforward way from the equations for the 
corresponding 2-point functions. Explicitly we use eqs. (13.42Q , Q3.50D , (|3.47|) , (|3.73|) and 
( |3.88 ) and obtain 



{k-D x )Kf,{k,X) 



(4.i; 



k ■ D x , J&(k,X) =2gNfk,k p F a (X)d»A{k) 



(4.2) 



k-D x ,J?Jk,X) 



= igk^ {^(X)t a Jf(k,X)-lfl(k,X)t a ^(X)} 

z^Nk^{^(X)lf a (k,X)-T/[ a (k,X)^(X)}, (4.3) 



k- D x , J^(k, X)] a =i 9 -Nk^ [i,{X) F(k, X) - %\k, X) i>(X)} 



(4.4) 



k ■ D x , J h \(k,X)] a = gN(d- 2) k^F^X) d v A(k), (4.5) 



(k-d x )B^k,X) = igk^ {xP(X))jL(k,X)-1/C(k,X)^(X)}, (4.6) 



(k ■ d x ) B 5 Jk,X) 



ig 



kfi {4>{X) l ,UKX)+1fi{k,X) l ^{X)}. (4.7) 



The equations ( |4. 1|) ( ^77|) are all first order partial differential equations and, to be spe- 
cific, we shall solve them in this section with vanishing initial conditions at X — > — oo. 
They describe the dynamics of the hard particles of the plasma in the presence of given 
background fields A^X), ip{X) and ip{X) which vanish adiabatically in the remote past. 

Only the real hard particles which exist in the plasma at thermal equilibrum do 
contribute to the induced sources. This is manifest in several ways in the previous equa- 
tions. First, in the r.h.s. of eqs. (|4.1|) - (|4.6|) , all possible vacuum contributions cancel. For 
instance, we can use the identity 



A(fc) + A(k) = Po (k)[N{k ) + n{k )} = 2ix5{k 2 )[N{e k ) + n{e k )}, 



(4.8) 



where e k = \k\, in order to see that the right hand side of eq. ( |4.1|) vanishes at zero 
temperature; the zero temperature contributions to the quark and gluon propagators 
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Sy^k) and Dy(k) are mutually compensating in their sum. Furthermore, the vacuum 
contribution to eqs. (|4.2|) and ( |4.5|) cancels out, as k • F • dk[(l/k)5(k ± e^.)] = 0. Thus 
only thermal fluctuations contribute to If. or J . Second, there is a factor (d— 2) (eventually 
implicit) in the r.h.s. of all the equations involving bosonic degrees of freedom. As already 
dicussed, this reflects the fact that only the transverse gluons effectively contribute to these 
densities. Third, the equations are independent of gauge fixing, which again reflects the 
decoupling of the spurious degrees of freedom. This is shown explicitly in Appendix B 
where the equations for the 2-point functions are studied for a covariant gauge with an 
arbitrary gauge fixing parameter A. 

All the previous equations are covariant under local gauge transformation of the 
mean fields A^, ip and ip. Therefore the induced densities determined from these equa- 
tions are gauge covariant quantities. Specifically, JjC transforms like ip as vector in the 
fundamental representation, while J' i A , jf y , and J A transform like F^ u as vectors in 
the adjoint representation, and B and B 5 are scalars. 

The similarity between the two equations for the vector and axial baryonic currents 
densities reflects the helicity conservation in the elementary interactions. The same sym- 
metry is, of course, valid for the corresponding color currents. The densities for the axial 
color currents satisfy equations which can be easily obtained by inserting 75 matrices in 
the vector equations ( |4.3|) and ( |4.4j ), in the same way as in going from eq. ( [4.6| ) to 

Quarks and gluons respond similarly to soft gauge and fermionic fields. This is 
apparent in the r.h.s. of eq. Q4.1|) , and also in the similarity between eqs. (|4.2| ) and ( |4.5|) . 
The different numerical factors in the right hand sides of these two equations simply 
count the relevant degrees of freedom. In fact, by adding together eqs. ([4.2|) and (|4.5|) , 
one obtains an equation for the total density J A = J A + J A 

k ■ D Xl J*(k, X)]=g k^k- F(X) ■ d (2N { A(k) + N(d - 2) A(fc)) . (4.9) 

This equation generalizes the linearized Vlasov equation to nonabelian plasmas. Note 
that previous attempts to derive such an equation were based on approximation schemes 
mixing leading and non leading contributions in g and therefore, from the present point of 
view, not entirely consistent |21, p2|. One can also combine eqs. (O) and (f4.4|) into a single 
equation for the total current density induced by the fermionic fields, = jf + , 

k-D x , J^{k,X)] =igk^t a {4>(X)t a ljC(k,X)-lfL(k,X)t a ij(X)}. (4.10) 



This equation is similar to the corresponding one in the abelian case [16|. In doing the 



sum of eqs. ( |4.3|) and fl4.4|) , the typical non-abelian effects cancel; these are contained for 
example in the second braces in eq. (fO) , and involve the 3-gluon vertex leading to gauge 
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field insertions on the hard gluon lines. This kind of cancellation was first noted by Taylor 



and Wong |26j in relation with the HTL's for amplitudes involving one pair of quarks and 
any number of soft gluons (albeit their proof is only explicit up to three external gluons). 

The right hand sides of the kinetic equations ( |4.1| )- (|4~7l ) contain a factor 5(k 2 ) show- 
ing that the mass shell condition for the hard particles is not affected by their interactions 
with the mean fields. This allows us, in the next subsection, to transform the equations 
for the current densities into kinetic equations for generalized distribution functions. 



4.2 Explicit solutions 

The solution of eq. (^4.1|) for has the following general structure: 

j£(fc, X) = 2ix5{k 2 ) {6(k )fi + (k, X) + 6{-k )4L-(-k, X)} . (4.11) 



This structure is characteristic to all the other densities that we compute below. Eq. ( |4.1| ) 
can be turned into equations for the functions JJL ± . It is easily verified that = = JJL, 
where satisfies 

(v ■ D x )0,X) = -i 9 -{d- 2)C f [N{e k ) +n(e k )]^(X), (4.12) 

— # 

with the initial condition fi.(k,X) — > for X — > — oo. Here j> = v^j^, with the velocity 
4- vector = (1, v) and v = kje k . This equation describes fluctuations where, under 
the action of a soft fermionic mean field, hard quarks are converted into hard gluons 
and vice versa. The quantity JjL(k,X) may be viewed as a generalized one body density 
matrix mixing fermionic and bosonic degrees of freedom. The covariant line derivative 
in the l.h.s. of eq. (|4.12|) is remininscent of the familiar drift term of elementary kinetic 



equations; the presence of a covariant derivative steems from the fact that JJL behaves as 
ip under gauge transformations. 



The solution to eq. ( 4.12 ) can be written by introducing the parallel transporter 



along a straight line joining x and y 

U(x,y) = Pexp{-z 5 / dz^A^z)}, (4.13) 

where 7 is the path parametrized by z^(s) = + (x M — y M )s with < s < 1, and P is 
the path-ordering operator. Thus, 

U(X,X-vu) = Pexp{-igu t dsv ■ A(X - vu(l - s))}. (4.14) 

Jo 
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In order to solve eq. ( 4.12|) , we take = A a ^t a in the fundamental representation. We 
have then 

Mk,X) = -i-(d-2)C f [N{e k ) +n{e k )]j) f°° duU(X, X - vu)^(X - vu). (4.15) 
2 Jo 

In order to verify that ( [4.15 ) is the correct solution to eq. ( 4.12|) we may use the following 
formula for the line-derivative of the parallel transporter 







(4.16) 



(vD x )U(X,Y) 

v A ' v ' ' Y=X-vu 

together with the identity (v ■ dx)f(X — vu) = — (d/du)f(X — vu). 

According to eq. (|3.8| ), the fermionic induced source is obtained by integrating the 



(4.17) 



density Jfi,(k,X) over k. By using eqs. ( [4. lip and ( |4.15| ), we get 

dVL 



7] md (X) 



, j) duU(X,X -vuWX -vu), 

4:71 JO 



with ojI = {d- 2)C f g 2 T 2 /16. 

The equation (|410|) for J* can be integrated in a similar fashion. One gets 

roo _ 

J* a (h, X) = igv'^ J o du U ab (X, X - v'u) 

{ip(X - v'u)t%(k, X - v'u) - Tft(k, X - v'u)t b ip(X - v'u)} , (4.18) 

where v' = k^/k = (1, ±v) for k = ±e k . U ab (X,X — v'u) is the parallel transporter in 
the adjoint representation, i.e., eq. ( |4 . 1 4] ) with — > = A a ^T a . The function Tfi has a 
structure similar to ( |4.11| ): 



Tfl(k,X) = 27T5(k 2 ){6(k )fi(k,X) + 6(-k )fi(-k,X)} 
where Uk,X) = ^{k,X)^°. 



(4.19) 



By integrating eq. ( |4.18| ) over k in order to get the induced current j^(X), we obtain, 
after simple manipulations, 



d 3 k 



(27T) 



3 ^ 



5n% a (k,X) -Snt a (k,X) 



(4.20) 



where 



5nt a {k,X) = ±i 



2e k Jo 



du U a b(X, X — vu) 



[i/j(X - vu)t b fi(k, X -vu)- A(/c, X - vu)t b ^j(X - vu)} 



(4.21) 
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are fluctuations in the quark and antiquark color densities induced by the fermionic mean 



fields. By replacing by its explicit expression ( 4.15 ), we can evaluate the integral over 
k = €k in ( [4.201 ) an d get 

/(ICy roo roo 

— v„ I du I du'U ab (X,X-vu) 



Jo 



i/j{X - vu)t h i}U(X -vu,X- v(u + u'))ip(X - v(u + u')) 

+^){X - v(u + u ))i)U{X - v(u + u'),X - vu)t b ip(X - vu)\. (4.22) 

We recognize in this expression the two terms which correspond to the two possible time 
orderings of the scattering processes between the hard particles and the soft fermionic 
fields. One can obtain a more compact expression by performing simple changes of vari- 
ables. Setting s = u + u' and t = u one transforms the first integrals over u, u' into 



oo roo 



dt / ds U ab (X, X - vtWX - vt)i;t b U(X -vt,X- vs)ip(X - vs). (4.23) 
o Jt 

In the second integrals, we set s = u and t — u + u', and obtain 

roo rt „ _ 

/ dt dsU ab (X,X -vt)ip(X -vt)4U{X -vt,X -vs)t b ip(X -vs). (4.24) 
Jo Jo 

By using the identity 

U ab (X,X -vt)t b = U(X -vt,X)t a U(X,X -vt), (4.25) 
together with the group property of the parallel transporter, 

U(X - vt, X - vs) = U(X - vt, X)U(X, X-vs), (4.26) 

we finally obtain 

/rjC) roo roo 

i)(X - vt)i;U{X - vt, X) t a U(X, X - vs)ip(X - vs). (4.27) 



We consider now the equations ( [4.6| ) and ( |4.7| ) for the baryonic currents. The 
solution to eq. ( |4.6|) is of the form 



B^k, X) = 2 {2n5(k 2 )) {6{k°)5n + {k, X) + 6(-k°)5n„(-k, X)} , (4.28) 

where Sn + (k, X) and 5n_(/c, X) are (color singlet) fluctuations in the quark and antiquark 
number densities. These satisfy 

(v ■ d x ) 5n±(k, X) = ±i^- \$(X)0, X) - X)^(X)] . (4.29) 



2e 



k 
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This equation is solved similarly to the previous one for J"^ (see eq. ( |4.21| )). Then we can 
calculate the baryonic current 

d 3 k 



and obtain 



d^k 
W) 



B,(k,X) 



4 M 



(27T) 



3 ^ 



5n+(k,X)-8n4k,X) 



(4.30) 



9U °J^Jo dt L ds 

ip(X - vt)jU{X -vt,X- vs)il)(X - vs). 



(4.31) 



This is, of course, similar to jff, eq. (|4 . 2 7| ) , but without the color matrices t a . In a complete 
analoguous way, we also obtain the axial baryonic current 

dVL 



/all f f 

ip(X - vt)^ U(X -vt,X- vs)ip(X - vs). 



(4.32) 



We turn now to the color currents induced by soft gauge fields. Eq. (|4.2|) for J A 
implies 



J^(k, X) = 2k ll N f t a (2n5(k 2 )) {6(k°)5n Aa (k, X) + 9(-k°)5n Aa (-k, X)} , (4.33) 

e 5n A = 5n Aa t a 
satisfy (with = \k\) 



where 5n± = 5n A . a t a are fluctuations in the quark color densities. These fluctuations 



v-D x ,6n A (k,X) = T gv-E{X) 



dn(e k ) 



de k 



(4.34) 



Here E(X) = E a (X)t a is the average chromoelectric field, E\ = F^°; possible chromo- 
magnetic contributions cancel because the equilibrium distribution n(e^) is isotropic. In 
the abelian case, this equation coincides with the linearized Vlasov equation. Here, the 
gauge field not only modifies the motion of the particle, through the color electric field 
that it generates, but it also induces a "precession" of the densities in color space, through 
the covariant derivative on the l.h.s. The solution of eq. Q4.34 ) is 



5n A {k,X) = T9 r duU{X,X - vu)v ■ E(X -vu)U{X -vu,X)^^- 
Jo dek 



(4.35) 



Similarly, eqs. (|4.5| ) or ( |3.93| ) describe oscillations in the gluon color densities 5N A = 
5N A t a . They satisfy an equation similar to (|4.34j) , viz. 

dN(e k ) 



v ■ D Xl 5N A (k,X) = -gv-E(X) 



de h 



(4.36) 



39 



The total induced current j A = j A + j A can be expressed in terms of these color density 
fluctuations. Thus 



j A (X) = gN { 



d 3 k 



(2*0 



3 



Sn A (k,X) -5n A (k,X) 



and 



jt(X)=g(d-2)N 



d 3 k 



(27T) 



VfM SN A (k,X). 



3 r 1 



(4.37) 



(4.38) 



Then, by using eq. (|4.35| ) for 5n±(k, X) and a similar expression for 8N (k, X), we obtain 

//•7Q poo 
—v» / rfw U(X, X - vu) F oj (X - vu)v j U(X - vu, X), (4.39) 
4:71 JO 

where = ((d — 2)N + N{)g 2 T 2 /18 is the plasma frequency. 

As a last remark, we note that the final expressions for the induced sources contain 
no explicit reference to the nature of the hard particles involved. This information is all 
contained in the explicit expressions of the two frequencies u and u p , which are obtained 
as simple momentum integrals involving the statistical occupation factors N(k) and n(k). 
These expressions may be easily generalized to the case of systems with a non vanishing 
chemical potential /i: 



ul = g 2 C f 



d 



16 



^ v = 9 l l-((d-2)N + N f )T 2 + ^ . (4.40) 



5 Mean fields dynamics 

In Sec. 3, we have derived consistent equations of motion for the 2-point functions which 
describe, at leading order in g, the response of the plasma to weak and slowly varying 
fields. By solving these equations in Sec. 4, we have expressed the induced sources, which 
enter the mean fields equations, in terms of the average fields themselves, thereby closing 
the system which determines the dynamics of the mean fields. 

The mean field equations, 

zp^X)=v(X)+V md (X), (5.1) 
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and 



[d\ f v ,{x) } a - g $(x)^ty(x) = j;(x)+jf da (x), (5.2) 

have a few remarkable properties, which are worth recalling. 

(a) As shown in Sec. 4.1, the induced sources are covariant under local gauge trans- 
formations of the mean fields. Specifically, rf nd transforms like if), while j transforms 
as F" . It follows that eqs. ( |5.1|) and (]5.2|) are gauge covariant if the external sources 
transform covariantly. This covariance is that of the Euler-Lagrange equations of motion 
deduced from the classical QCD action, (i.e., the tree level action, eq. (|2.4|), but without 
the ghosts and the gauge fixing terms). It is noteworthy to recover gauge covariance at 
the level of eqs. ( |5.1| )-( |5\2j ), where corrections to the tree level dynamics are taken into 
account through the induced sources. 

(b) The derivatives of the mean fields in the l.h.s. of eqs. (|5.1|) and fl5.2|) are 
of the same order of magnitude as the induced sources in the r.h.s. That is, r] md ~ 
gTip ~ dxip and j A ~ gTF ~ dxF, as can be readily verified on the explicit expressions 
given in eqs. ( |4.17| ) and ( [4.391) . Thus, the propagation of soft modes is nonperturbatively 
renormalized by the medium, reflecting the collective character of the long wavelength 
excitations. 

(c) For average fields whith maximum strength, i.e., F ~ gT 2 and ipip ~ 9T 3 , all 



the terms appearing in eqs. ( |5.1| )-( |5T2| ) are of the same order of magnitude. The equations 
are then nonlinear, due to the covariant derivatives in their l.h.s., as well as those which 
enter the structure of the induced sources. This is the only source of nonlinearity, showing 
that the limits on the fields strength we are using are equivalent to neglecting all genuine 
non linear dynamical effects of the fields. 

(d) Eqs. ( |5.1D and ( |5.2|) are nonlocal, which is a manifestation of the fact that the 
ultrarelativistic QCD plasma is a dispersive polarizable medium. 

For vanishing external sources, the mean field equations describe the normal modes 
of the plasma, or equivalently its soft quasiparticles. These equations have been solved 
only in the case where the gauge fields are weak, that is in the linear, abelian, regime. 
Then the oscillations of the fermion and gauge fields decouple, and explicit dispersion 
relations can be obtained 0, [15, |34|. In general, quark and gluon modes mix, and 



the non linearity of the equations make them difficult to solve. Leaving aside this difficult 
problem, we note that the mean field equations summarize a lot of information about the 
soft quasiparticles, such as their effective propagators, and the effective vertices which 
describe their mutual interactions. We show in the next subsections how these can be 
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obtained from the induced sources, and derive the simple Ward identities that they satisfy. 



5.1 Hard thermal loops from induced sources 

As pointed out in Sec. 2. 2, the knowledge of the induced sources as functionals of the fields 
allows us to derive, by functional differentiation, the one particle irreducible amplitudes 
for the corresponding fields. We use this procedure now, in order to construct self-energies 
and proper vertices for soft quasiparticles valid at leading order in g. In doing this, we use 
boundary conditions such that the mean fields vanish at Xq — ► — oo, as in the previous 
section. The resulting amplitudes will be causal with respect to the time argument in the 
induced source. 

In order to proceed in a systematic way, it is convenient to introduce retarded and 
advanced Green's functions for the covariant line derivative v-Dx, denoted as G ret (X, Y; v) 
and G a d v (X, Y; v), respectively. Thus, G ret (X, Y; v) is defined by 

i (v ■ D x ) G ret (X, Y; v) = 6 W (X - Y), G ret (X, Y; v) = for X < Y , (5.3) 

and G a dv{X, Y;v) satisfies a similar equation, but G a d v (X, Y; v) = for X > Y . The 
solution to eq. (|5.3|) may be written as 



G ret (X,Y;v) = -i6(X° -Y°)5 (3) (X -Y - v(X° -Y ))U(X,Y) 

= -i / du 5^ (X — Y — vu) U (X, Y) , (5.4) 
Jo 

where U(X,Y) is the parallel transporter ( (4.13j ). The use of this Green's function is 
appropriate for the study of the response of the plasma to perturbations which vanish 
as Xq — > — oo. One could also consider perturbations which vanish adiabatically when 
X — > oo. Then, the induced currents would be expressed in terms of the advanced 
Green's function, 

roo 

G ad JX,Y;v) = i / du5 {4) (X-Y + vu)U(X,Y), (5.5) 
Jo 

which is related to G ret by G a d v (X, Y\ v) = G* ret (Y, X; v), where the Hermitian conjugaison 
refers to color indices. For ^4^ = 0, the Fourier components of G re t an G a dv are simply 
G ret (P) = (v-P + tr ] )- 1 = G a * dv (P), with n -+ 0+. 



According to equations like ( p. 271 ) or (|2.30|) , r] md acts as a generating functional for 



all the proper vertices having at least one pair of external fermion lines and any number 
iV of bosonic lines. In the present approximation, r] md is given by eq. (|4.17|) , which can 
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also be written as 

v ind {X) = J d A Y5J: A {X,Y)tfj{Y). (5.6) 

Since r] md is linear in ip, the only amplitudes which do not vanish at the level of the present 
approximation involve in fact no more than one pair of external fermions. The kernel 5T* A 
may be identified with the self energy of a soft fermion propagating in a background gauge 
field (see eq. ( |2.27| )). In order to get its explicit form, we rewrite eq. (|4.17|) for r\ md as 

dVL 



r jind ! 



\X) = u 2 J — ^Jd 4 YG ret (X,Y;v)^(Y). (5.7) 
Then, obviously, 

6Z A (X,Y) = u 2 J^^G ret (X,Y;v), (5.8) 



which is similar to the corresponding abelian expression [IB]. The only manifestation of 
the underlying gauge group is the presence of the color factor Cf in the expression of a>o, 
together with the non-trivial color structure of G re %. Eq. ( |5.4j ) implies <5S A (X, Y) — for 
X° < Y°, reflecting the retarded boundary conditions in eq. ( |5.7| ) for rf nd (X). For A = 0, 



eq. (|5.8| ) gives the well-known HTL expression for the soft quark (retarded) self energy 



mf) = ^ I — — i — . (5.9) 



Note that, because P ~ gT, the correction (|5.9|) to the quark inverse propagator is of the 
same order of magnitude as the tree-level amplitude (Sq 1 (P) ~ If ~ gT). It can be used 
to construct the equilibrium propagator for a soft quark, correct at leading order in g: 
*S- 1 (P) = -f + 5Z(P) §. 

Because eq. ( |5.8| ) is valid for arbitrary (soft) background gauge fields, we can use 
it to derive - by successive functional differentiations - expressions for all amplitudes 
involving one fermion pair and any number of soft gluons. To be specific, we define the 
correction to the amplitude between a quark pair and N soft gluons by 

g N 6r%;;Z(X 1 ,...,X N ;Y 1 ,Y 2 ) = — — - ~ — — - (5E A {Y U Y 2 )) . (5.10) 



5A!£(X N )...5A£l(X h 



In doing these differentiations, we use identities like 

# G ret (Xx, X 2 ; v ) . . . / C ii\ 
SA ^ Y ) = 9v^G ret ( y X 1 ,Y;v)t G ret (Y, X 2 ; v), (5.11) 

and similarly for G a d v , which follow from eqs. (|5.4j) and (|4.13|) . The normalization we 
choose for the amplitudes ( p.!0[ ) is such that 5T depends on g only through uJq. For 
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all these amplitudes, Y® is the largest time, while Y 2 ° is the smallest one. The relative 
chronological ordering of the N gluon lines is arbitrary, and, in fact, the amplitudes 
are totally symmetric under their permutations. Up to minor changes due to the color 
algebra, all the amplitudes obtained in this way coincide with the corresponding abelian 
amplitudes |16| . 



For the sake of illustration, we give now the explicit expressions for the amplitudes 
involving a quark pair and one or two gluons. The first case corresponds to the quark- 
gluon vertex correction in the presence of a background gauge field (see eq. (|2.30|) ) 

5V A ;{X-Y U Y 2 ) = w o y j^v ll v v G ret (Y 1 ,X;v)t a G ret (X,Y 2 ;v), (5.12) 



where the time arguments satisfy Y± > X° > Y 2 °. For A = 0, eq. ( |5.12|) gives the HTL 



correction to the quark-gluon vertex p|, with appropriate boundary conditions. Defining 
the Fourier transform of <5r^ by 

(2vr) 4 c^(P + K x + K 2 ) SY^P- K u K 2 ) = 

J d*Xd A Y t d A Y 2 e i(P^+Ki-yi+K 2 -Y 2 ) 6T a^ x . ( 5 _ 13) 



we get 



(5.14) 



As the external momenta are all of order gT, g5T^ is of the same order as the bare 
vertex gj^. Thus, the effective quark-gluon vertex at leading order in g is gt a *r^, where 

*r M = TM + av §. 

The second example is that of the vertex between a quark pair and two gluons, 
g 2 *T a ^ u . This vertex does not exist at tree level, and in leading order it arises entirely 
from the hard thermal loop, *T°^, = ST"^,, with 
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STf u (X 1 ,X 2 ;Y 1 ,Y 2 ) 



uj y i — Vfl v v v p 



G ret {Y x , X i; v) t a GretiXt, X 2 - v) t b G ret (X 2 , Y 2 - v) 
+G ret (Y u X 2 ; v) t h G ret (X 2 , X x - v) t a G ret (X u Y 2 ; v)}. 



(5.15) 



For A = 0, defining the Fourier transform as in eq. (|5.13|) , we obtain 

dVL 



47r (v ■ K\ + i7])(v ■ K 2 — irj) 

t a t b t b t a 



+ 



(5.16) 



[ v ■ (K x + Pi) + irj v {K x + P 2 )+iri}' 
which, when summed over the color indices, coincides with the expressions given in refs. ||, 



Alternatively, we can derive the amplitudes ( 5.1 0|) from the expression ( 4.271 ) for 
the induced current j^(X) = jf(X) + j^(X), but, this time, the boundary conditions 
will be different, as the time argument of j^(X) is now the largest one. We first use the 
definitions (|5.4| ) and (|5.5|) to rewrite j^(X), eq. (|4.27|) , as 



j*(X) = gt^l f d ^-v, J d^Y 2 



^(Yi) i G adv (Y u X- v) t a G ret (X, Y 2 ; v) ^(Y 2 ). (5.17) 
The resulting correction to the quark-gluon vertex is (recall the second equality in eq. ( |2.30| )) 

5Tj a (X;Yi,Y 2 ) = J ^ v^v u G adv (Y 1} X; v) t a G ret (X, Y 2 ; v), (5.18) 

where now the time arguments satisfy X° > max(Y®,Y 2 °), the chronological order of Yi 
and Y 2 being arbitrary (compare, in this respect, with eq. (|5.12|) above). For A = 0, we 
have 



STliP-K^K^, 



j.a, ,2„,f 

-t o; 7 



dQ 



An (v • K\ — ir/)(v ■ K 2 — irj) 



(5.19) 



which differs from ( |5.14|) solely by the irfs in the denominators reflecting the respective 
boundary conditions. If we further differentiate the 5T^ a (Xi; Y 1; Y 2 ) of eq. ( |5.18|) with 
respect to A^, we generate amplitudes of the type ( |5.1(J| ), in which Xf is the largest time. 
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All the amplitudes involving only soft gluons may be derived from the induced 
current j A . In the present approximation, j A is given by eq. (|4.39|) , which, after some 
algebraic manipulations, can be cast into the form 

r r/O r 

j Aa (X) = -3^A a (X)g , + 3i^j —v,j d A YGf et (X,Y)(v.A b (Y)), (5.20) 



where G re t is the retarded Green's function ( |5.4|) written in the adjoint representation 
(i.e., with U —>■ U), while A = d$A. By successive differentiations with respect to A 
in j A (X), eq. (|5.20|) , we derive equilibrium amplitudes with any number of soft gluon 
external lines, which are causal with respect to Xq. For example, by using eq. ( |2.26[ ) one 



SK(P) = 3^ -§y u + P° I \ . (5.21) 



gets, after a simple calculation, the well-known expression of the (retarded) soft gluon 
polarization tensor: 

Air v ■ P + if] 

As P ~ gT, this is of the same order as the tree-level inverse propagator, Dq X (P) ~ 
P 2 ~ g 2 T 2 . Thus, in leading order, the effective propagator for soft gluons is *D~ 1 (P) = 
Dq^ v (P) + 5H^ V (P) m, 0. Note that the constant term in eq. ( |5.21|) stems from the first 



term (linear in A ) of the expression eq. ( |5.20p for j A ; this term does not contribute to 



amplitudes with N > 3 external lines. As a further example, we write the expression for 
the three-gluon vertex correction which follows by a second differentiation with respect 
to A (as in eq. ( |2.31| )). Using the same definition as in ( [5.13D for the Fourier transform, 
we obtain from eq. Q5.2C ) 



SFf^P^Ps) = rf abc 3u 2 p J 



abcn,2 V ^uV p 



pO pO 



Air v ■ Pi + ir] 

■ \ , (5-22) 
I v - "3 — irj v ■ P2 — ir] J 

where the imaginary parts in the denominators correspond to the time orderings A° > 
X® > A3 for the first term inside the parantheses, and, respectively, A° > A3 > A° for 
the second term. We can rewrite this more symmetrically as 5T ab ^ p = if abc 5T M1/p with |3j 

2 r cM f P° - P 2 ° 



ST IMl/p (P l ,P 2 ,P 3 ) = u 2 J —v^Vp 



(v ■ Pi + irj)(v ■ P 2 - ir]) 



pO _ pO pO _ pO ~| 

+ ^ j3 + n n 1 / 5 _ 23 n 

(v ■ P 2 - ir/) (v ■ P 3 - irj) (v ■ P 3 — ir]) (v ■ P x + irj) J ' 

Generally, all the amplitudes with soft external gluons which are derived in this 
way coincide with the corresponding HTL's of the diagrammatic approach |3|, 0, with 
boundary conditions such that A° is the largest time, while all the other external lines 
enter symmetrically. 
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5.2 Conservation laws and Ward identities 



The induced currents satisfy covariant conservation laws related to the symmetries of the 
QCD lagrangian. These can be easily derived from the equations of motion for the corre- 
sponding densities. These conservation laws imply, in particular, simple Ward identities 
relating the hard thermal loops. 

By acting on eq. ( f4.10|) with the covariant derivative D x , we get 
D x , [k-D x ,J*(k,X)\] = [k ■ D x , igt a (4j{X)t%{k,X) -Tfi{k,X)t a ^j{X))] .(5.24) 
In the left hand side, we use the Jacobi identity to write 



u Xi 



k ■ D x , J L 



7 <-> fj, 



k-D 



x ■ 



J;, F^K 



(5.25) 



where the definition of the field tensor igF^ = [D^, D u ] has been used. Since i 



is 



proportional to k^, as shown by eq. Q4.10D , the last term cancels in the right hand side of 
eq. ( ggg ). Then eq. ( gg§ leads to 



D x , Jt{k, X)] = ig t a ($(X) t a UK X) - nk, X) t a ^(X)) . 
After integrating over k, this becomes a continuity equation for 

D\ j*{X)] = igt a (${X)t a rf nd {X) - ff nd (X)t a ^(X)) . 



(5.26) 



(5.27) 



Furthermore, the equations of motion for the fermionic fields ip and ip, that is, eq. ( |5.1| ) 
and its Hermitian conjugate, may be used to show that 



igt a [ifjt a r] - rjt a ip 



ig t a (i)t a 7] ind - ff nd t a 



(5.28) 



By adding eqs. (|5.27|) and (|5.28|) , we finally obtain the continuity equation satisfied by 
the total current gip'j fl t a ip + associated to the fermionic mean fields: 

D», gty^^ty + jf\ = -igt a {^Jt a r] - r/t>) . (5.29) 

Thus the current gipj^ip + is conserved in the absence of external fermionic sources. 

Similar conservation laws hold for the axial color current and for the baryonic cur- 
rents ( [OflD and flP2p . 

Repeating the same steps as above, but starting with the Vlasov equation ([4.9|) , we 
obtain the covariant conservation law for the induced current j A : 



o. 



(5.30) 
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Finally, by using the conservation laws for the induced currents, eqs. ( |5.29| )-( [5.30| ), 
together with the Jacobi identity [D u , F v/i ]\ = 0, we see that, in order for the mean 
fields equations ( |5.1| ) and ( |5.2| ) to be consistent, the external sources must satisfy 

[D^j lx ]=igt a (i>t a i 1 -f\ty). (5.31) 

By differentiating the equations ( |5.30|) and ( |5.27|) with respect to the fields, one 
obtains relations between HTL refered as "QED-like" (or "tree- level like" ) Ward identities 
in Refs. |3|, [J. In particular, the transversality of the polarization tensor follows from 
eq. (|5.30|) , after one differentiation with respect to A^. 



m^P" = p^sn^ = o. (5.32) 

Of course, this property can be easily verified on the explicit expression for SU^ given in 
eq. (|5.21|) . If we differentiate once again eq. (|5.30|) , we derive the simple identity 



P?8T^ X {P X , P 2 , P 3 ) = 5U uX (P 3 ) - 5U uX (P 2 ), (5.33) 
which also follows from eqs. (|5.22|) and (|5.21 ) above. 



For amplitudes involving a quark pair, Ward identities result from differentiating 
eq. ( |5.271 ). For instance, by differentiating eq. ( |5.27f) with respect to both if) and if), and 
using eqs. ( |2.27|) and (|2.30|) , we recognize the QED-like Ward identity 

P^r M (P; K u K 2 ) = 5Z(Ki) - SE(Kx + P), (5.34) 



which can also be deduced from the explicit expressions eqs. ( p.!4| ), ( |5.19| ) and (|5.9|). As a 



final example, consider the following identity, which results from eq. ( |5.27 ) after an extra 
differentiation with respect to A\, 

Pr5Tf u (P 1 ,P 2 ;K 1 ,K 2 ) = if^TliP. + P^K^K,) 

+ 5T b l/ (P 2 ;K u K 2 + P 1 )t a - t a 5T b u {P 2] K 1 + P l ,K 2 ), (5.35) 

and can also be obtained directly from eqs. (|5.14 ) and ( 5.16|) . 



5.3 Effective action for soft fields 

Using the results of Sec. 5.1, one can rewrite the equations of motion for the average fields 



as 

dVL 



— f J d A YG ret (X,Y;v)i)(Y) = V (X), (5.36) 
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and 

[ D", F Ufl (X) f - 3 zul J — v, J d A Y G? et (X, Y) v ■ E b (Y) 

= j;(X) +^(X)W(X) + j; a (X), (5.37) 

where jt{X) is the current induced by fermionic mean fields, given in eq. (pTj) . In these 
equations, the fields are assumed to vanish adiabatically for Xq — > — oo, which correspond 
to the plasma being initially in equilibrium. Of course, similar equations may be written to 
describe perturbations which vanish adiabatically in the remote future; then, the induced 
currents would be expressed in terms of the advanced Green's function. 

In order to make contact with previous works, we shall now restrict the fields ip{X), 
4>(X) and A^X) to belong to a space TZ in which the operators v ■ Dx, v ■ Dx and 
(v ■ Dx) 2 never vanish. Specifying the boundary conditions becomes then unnecessary, 
and it is possible to deduce the equations of motion ( |5.36| ) and ( |5.37| ) by applying the 



minimum action principle to an effective action which can be written in a simple compact 
form. The space 1Z is spaned by the fields which vanish exponentially for \Xq\ —* oo and 
satisfy the conditions 

/oo 
du U(X, X - vu)tj)(X - vu) = 0, (5.38) 
-oo 

/oo 
duU ah {X,X-vu)Fl v {X-vu) = 0, (5.39) 
-oo 

/oo 
duuU ab {X,X -vu)F^(X -vu) = 0. (5.40) 
-oo 

These conditions guarantee that the fields have no components in the kernel of the covari- 
ant line derivatives. The left hand sides of eqs. (|5.38|) - (|5.40|) represent such components. 
For instance, 



(vD x ) du U{X, X - vu)ip{X - vu) = 0, (5.41) 

J — oo 

as one can easily verify (see e.g. eq. ( |4.16|) ). Physically, the conditions above imply that 
the fields do not suffer Landau damping. This is seen most simply in case of weak fields, 



where we can ignore the parallel transporter. Then the Fourier transforms of eq. ( |5.3£ 
and (|5.39| ) result in 



ifj(P° = v . p, p) = 0, and F* U (P° = v ■ P, P) = 0. (5.42) 
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Within TZ, the inverse operators \jv ■ D, 1/v ■ D and l/(v ■ D) 2 are nonambiguous and 
represent indifferently the retarded or advanced Green's functions, which are equivalent 
on TZ. 

We write the effective action as S e ff = So + Si n d, where So is the classical QCD 
action, (as given in eq. (|2.4j ), but without the ghosts and the gauge fixing term), while 
Si n d contains the effects of the interactions between the soft fields and the hard particles. 
In order that S e ff generates the mean field equations of motion ( |5.36| ) and (|5.37| ), one 
requires 

5S ind /5$(X) = rf nd {X), 5S ind /5A»(X) = f a nda (X). (5.43) 

As proven in Appendix C, these conditions are satisfied by Si n d = Sf + Sb, with 

S t = ^/^/rf 4 x|d 4 F^(X)(A|^i^|r)^(r), (5.44) 



and 



F^(X)(X\-^-\Y)F-\Y) 
[v ■ DY 



(5.45) 



where the trace acts on color indices only. The action St accounts for all the polarization 
effects induced by soft fermionic mean fields, that is 

8St/8i>{X) = rT*(A), 8St/6A£(X) = j% a (X), (5.46) 

while Sb satisfies 

5Sb/5A»(X) = rf a {X). (5.47) 



The effective action ( 5.44f )- (|5.45|) coincides with the generating functional for HTL's first 



derived in [^, ^3 on the basis of gauge invariance. One sees from our analysis that 
S e ff may be given the physical interpretation of the classical action which describes long 
wavelength excitations of the hot quark-gluon plasma, at leading order in the coupling 
g, but only as long as Landau damping is inoperative. One should emphasize that, 
indeed, the functional derivatives of Si n d generate properly the induced currents only on 
TZ (see the discussion in Appendix 3). In cases where the boundary conditions matter, 
we have to recourse to the kinetic equations written in Sec. 4.1 to solve for the induced 
sources appropriate to the given physical conditions (as done e.g. in Sec. 4. 2 for retarded 
conditions.) 
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6 Conclusions 



We have presented a consistent description of soft collective excitations of quark-gluon 
plasmas. This description relies on the hierarchy of Dyson- Schwinger equations, which in 
the leading order that we have considered, truncate at the level of 2-point functions. Our 
resulting equations take then the form of a coupled set of mean field and kinetic equations. 
The mean fields describe the soft excitations and have many features of classical degrees 
of freedom. The kinetic equations describe the dynamics of the hard plasma particles, 
which respond collectively to the action of the mean fields. This separation between mean 
fields and kinetic equations reflects the distinct roles played by the two types, hard and 
soft, degrees of freedom of the plasma, and also the fact that the dominant interactions 
at leading order are those between the hard particles and the mean fields. 

As hard particles behave in a remarkably similar fashion, irrespective of whether 
they are quarks or gluons, they provide sources not only for average gauge fields, but 
also for average fermionic fields. And indeed the quark-gluon plasma sustains collective 
excitations with fermionic quantum numbers. That is, there exists excitations in which 
the soft fermion degrees of freedom oscillate over a long time scale. These oscillations 
are maintained by the hard particles coherently changing their quantum numbers from 



bosonic to fermionic and vice versa, and this over long space/time scales [34 . 

The excitations that we have studied are related to oscillations in the conserved 
currents associated with the QCD Lagrangian. There are thus oscillations in the baryon 
currents which result from the action of a fermionic mean field. There are also oscillations 
of color currents induced either by a fermionic field or by a gauge field. All these oscil- 
lations can be described in terms of fluctuations of density matrices which are generally 
non diagonal either in color space or in a space containing fermionic and bosonic degrees 
of freedom. The diagonal matrix elements, when they are non vanishing may be given 
the interpretation of classical distribution functions for particles of a given color. 

The simple physical picture that we have just outlined is of course fully justified 
by the long technical developments presented in this paper, whose main features we now 
summarize. 

The cornerstone of the whole approach is that collective modes develop at a partic- 
ular energy scale, gT. The emergence of this scale is seen most clearly on the mean field 
equations which, in leading order, become universal equations, once length are measured 
in units of 1/gT and energies in units of gT. It is the existence of this particular scale 
which allows us to organize our analysis of the Dyson-Schwinger equations, by controlling 
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consistently at the same time the strength of the coupling, the wavelength and the am- 
plitude of the collective oscillations, to eventually arrive at a consistent set of equations 
of motion which incorporate all leading order effects. 

One outcome of our work is a set of consistent and gauge covariant kinetic equations 
for the plasma particles moving in background fields. As another by-product, we recover 
in a simple way all hard thermal loops which have been analyzed first in a diagrammatic 
approach. These are simply obtained from compact expressions for the induced sources, 
which may be regarded as generating functionals for the HTL's. Thus our approach 
provide a new and independent derivation of such generating functionals. It has also the 
advantage of clearly exhibiting the role of boundary conditions. In particular our induced 
currents are easily obtained with the retarded boundary conditions appropriate to the 
calculation of the causal response functions. The standard form of the HTL generating 
functional is recovered by restricting the fields to obey special conditions. The particular 
Ward identities which relate HTL's follow simply from conservation laws obeyed by the 
induced currents. Finally, the mean field equations provide expressions for the resummed 
soft propagators and vertices. 
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Appendix A 



We prove here that the propagation of soft mean fields with ghost quantum numbers 
does not give rise to polarization effects at leading order in g. That is, in this order, the 
induced sources can be ignored in the r.h.s. of eqs. (|2.20| ) and (|2.21|) for the ghost mean 
fields. (In the diagrammatic language, this is equivalent to say that there are no "hard 
thermal loops" with ghost external lines [0].) It follows that the ghost fields vanish in 
the absence of external sources and they can be omitted in the equations of motion for 
the 2-point functions (see Section 3). For brevity, we present the explicit proof only in 
Feynman gauge and for vanishing gauge fields. The general case can be treated along the 
lines described at length in Appendix B. 

The induced sources C™ d and C™ d are related to the diagonal matrix elements of 
the ghost-gluon off-equilibrum propagators, O ab and 0^ b , by eqs. ( |2.44| ) and ( |2.45| ), which 
imply 



and 



C mda (X) = gf 



abc 



dxO u 



bc<. 



s,X) 



s=0 



gf 



abc 



d A k 



d x o b ;(k,x) } 



(A.l) 



C mda (X) = gf 



abc 



s,X) 



s=0 



-igf 



abc 



d 4 k 



k"Ol c (k,X). 



(A.2) 



Note that in going from eq. (|2.45| ) to eq. (|A.2|) we have neglected the soft derivative dx 
in comparaison with d s , whereas ( |A.1|) is an exact rewriting of eq. fl2.44j ). 

The equations for O ab and are obtained by differentiating the mean fields equa- 
tions with respect to the external sources, as explained in Sec. 3. 3. First, we differentiate 
the gauge field equation (|2.18|) with respect to C b and C b , and get, respectively, 



££(y)0£(s,y) 



-gf 



acd 



die 



be i 



*,v)l <C d (y)>, 



(A.3) 



and 



y) = gf acd G cb (x,y)d fl (C d (x)). 



(A.4) 



Then, we differentiate eqs. ( |2.20| ) and Q2.2ip for the ghost fields with respect j^, and 
obtain 



di O ab (x, y) - gf acd (A c (x)) d»O db (x } y) 



gr d d» DZ(x,y){C(x)) 



/IV 



(A.5) 
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and 



d 2 O h v a (x, y) - gf acd (Afa)) 6%&«{x, y) 



gf acd D b u l(x,y)d»(( d (y)) (A.6) 



Let us consider first C™ d . This will be obtained after solving eqs. ([A. 3 ) and ( |A.5| ) for 



O ab . As mentioned earlier, we choose (A*) = 0, in order to simplify the discussion. Then, 



after doing approximations similar to those described at length in Section 3, we obtain 
for O ab the following two equations: 

d 2 O ab <(x,y) = 9 r bc (dyA(x-y))C(y), (A.7) 

and 

dl O ab< (x, y) = -gf abc dl [A(x - y)( c (x)} . (A.8) 

The ghost mean fields are simply denoted here by ( a and ( a . These equations are further 
transformed into an equation for the Wigner function O ab (k,X) 

(k-d x )O ab (k,X) =if /^M^CXqjA). (A.9) 

This equation includes consistently all the terms at leading order in g. It shows that 
0(k,X) ~ g((/T)A. One deduces then from eq. (|A~T| ) that C ind ~ g 3 T 2 (, which is g 
times smaller than the l.h.s. of eq. ( |2.20| ) for d 2 ( ~ g 2 T 2 (. In fact, C md is even more 
suppressed than this estimate suggests. 

The corresponding proof for the anti-ghost mean field - i.e., for C md and 0^ b - is 
slightly different, owing to the disymmetry with which the two kinds of ghost fields enter 
the QCD action ( |2.4| ). Now it is the abnormal propagator 0^ b itself which vanishes at 
leading order. This can be easily verified by isolating the dominant terms in eqs. ( A.4j) 
and for Of. 



Appendix B 

In the derivation of the equations of motion in the main text, we have used sys- 
tematically the Feynman gauge (A = 1) in order to simplify the algebra. We prove now 
that this particular choice does not affect the final results: the equations that we have 
obtained for the mean fields and their induced sources preserve their form in any covariant 
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gauge. As a consequence, the HTL's are independent of the gauge fixing parameter A for 
arbitrary (soft) external momenta. 

The gauge fixing parameter A enters our calculations through the gluonic prop- 
agator, that is, either through the differential operator C? v (x), eq. ( |3.18|) , or through 



Duuvik) = 5 p /ll/ (k)N(ko). Here p^k) is the spectral density in a general covariant 
gauge [H 

pAk)=P F ,Ak) + (l-\- 1 )pl{k), (B.l) 

with p^ v {k) = —g^uPo(k) the spectral density in Feynman gauge (po{k) = 2Tre(k )5(k 2 )), 
while Pp U (k) = —27ie(ko)kfj,k u S'(k 2 ), S'(k 2 ) denoting the derivative of 8{k 2 ) with respect 
to k 2 . We also define A'(k) = 2ire(k )5'(k 2 )N(k ), by analogy with A(k) = p (k)N(k ). 
By using k 2 5'{k 2 ) = —S(k 2 ), we verify that 

k 2 A'{k) = —A(k), (B.2) 

and consequently the gluon propagator can also be written as 

D^M = W»k 2 - (1 - X-^kMA'ik). (B.3) 

We shall need later the following identity: 

k»D< tlu (k) = -\- l k u A(k). (B.4) 

After these preparations, let us formulate the proposition to be proved below: in a 
general covariant gauge parametrized by X, the k-space densities for the induced sources, 
i.e. fyl(k,X), jf(k,X), J^(k,X), J~ { A (k,X) and J^(k, X), are independent of X when 
consistently computed at leading order in g. 

Note first that, since the gluon propagator is not involved in the equations for 
SS A , J A is trivially gauge independent. Furthermore, by examining the derivation of the 
equations satisfied by jf and jTjf, one sees that these may depend on A only through tfL 
and TfL. Thus, we only need to prove that )f and J A are independent of A. 

We start with If,. By following the same steps as in Sec. 3. 3, we derive from eq. Q33 
the equation which generalizes eq. (|3.33|) to arbitrary A 



k 2 + ik ■ (d x + 2igA)] « -gA(k) [jfy, - (1 - A" 1 )^] t a ip, (B.5) 

where we have preserved only the terms of relevant order of magnitude. Similarly, starting 
from eq. ( |3.28|) we get 



(k 2 - ik ■ d x )KZ ~ tgf abc T^ pX A^IC u b + i(l - \)hJC a « gA{k)^,ty. (B.6) 
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We have introduced here the function JC a (k,X), the Wigner transform of 

K a {x,y) = d v y K< a {x,y). (B.7) 
We now expand K, = /C (0) + /C (1) + with /C (1) ~ #/C (0) , etc. At leading order, we 



recover from the previous equation the consistency conditions ( |3.35|) and ( |3.39|) , that is, 
k 2 K,^ = and k^K,^ = 0. At next to leading order, we subtract the two equations and 
obtain the generalization of eq. ( |3.37| ) to arbitrary A 



k • (dx + igA) Kl + 9 -r bc T, up ^ c k x K b 



n rou,r^a 



.9 



i^-[A(k) + A(k)]^t^ - 



■9. 



A- 1 )A(A;)V>, 



(B.8) 



where /C = /O ). Because k^JC^ = 0, the leading order contribution to JC a (which enters 
eq. O) is ~ gTK<®: 



ic a {k,x) = ik»jcp a (k,x) + l -d x K^ a (k,x). 



(B.9) 



Thus eq. ( B.8|) is consistent, i.e., all terms are of the same order of mag nitude, ~ gT 2 JC(°\ 
The unknown function KS l \ as well as the A dependence disappear in the equation for 
Ifi 1 = J^fcffl. These cancellations result from the identity 



$K a {k,X) = -i^A(k)t a ij(X), 
A 



(B.10) 



which is proved by first differentiating eq. ( |3.29|) with respect to dy, and then using the 
identity (pip , to get 



s K a (x, y) = -i^n a ^)d^A{x - y), 



(B.11) 



which, in leading order in g, gives precisely eq. (|B.10|) . Thus, we recover for )f a the 
Feynman gauge equation (|3.40|) . 

Consider now the gluonic current j£. As discussed in Sec. 3. 7, = + Jg + Jo,- 
Since, by construction, J G is independent of A, we need to prove that a similar property 
holds for the sum + Jq. 



From eq. ( 3.15 



j^(X) = 9 2 f;5^(I)^(0) = Ng 2 [g, u (D<r p (0) - (D<), u (0)\ K(*)- ( B - 12 ) 

By using eq. ( |B.3| ), this can be written as a integral in /c-space of Jq = j£ + (1 — A _1 )j7 A , 
where 



q-F a 



-gN(d-l)A(k)A a JX), 



(B.13) 



56 



is the Feynman gauge expression, eq. ( |3.17| ), while 

J x ^(k,X) = -gN(g^k 2 - k^K)A'(k)A:(X) (B.14) 

is a supplementary term for A 7^ 1. 

As indicated in eq. ( |3.76| ), is determined by the quantities SV pu k u and bV> a v v . 
(To simplify the writing, we omit, here and below, the A-superscripts.) These quantities 
will be determined by following the same steps as in Sec. 3. 7. We start from eqs. (|3.58|) 
and ( |3.59|) , written for arbitrary A and ip = ip = 0. By using eq. (|3.18|) for Cw v , we write 
eq. (|3.58|) as 

(d 2 X " (1 " A)^) Df v <(x, y) = -gf adc [r ^ pX A>(x)d x x - T, pXa (d x x A p c (x))] D£<{x, y) 

+ g 2 t<^ pX Al{x)A*{x)D£<{x,y). (B.15) 

We shall refer to this equation as (I). We have a similar equation deduced from (|3.59| ) and 
to which we refer below as (II). It is obtained from (I) by exchanging, respectively, x and 
y, a and b, and fi and v. We separate the induced part of D by setting D < = D$ + 5D, 
with 3D ~ gD for A ~ T. Note that D$ does not contribute to the l.h.s. of (I) and 
(II), because it satisfies Dq 1 Dq = 0; (Dq 1 is the free inverse propagator, eq. (|2.28|) ). We 



expand SD = 5D^ + 5D^ + where <5L> (1) ~ gSD^°\ etc... At leading order, (I) and 
(II) lead to the consistency condition which generalizes eq. ( |3.61|) for A 7^ 1 

k 2 6V a ^ - (1 - X)k^5V a pi/ = gT, upX AP c (X)k x D<:(k), (B.16) 

where we have set (see eq. (|3.74| )) 

5V; u (k,X) = (i/N)r bc 8V$ bc (k,X). (B.17) 

At next-to-leading order, we encounter terms involving the unknown function 5D^\ 
For instance, (k ■ d x )6D$ ~ k 2 5D$ ~ k p k p 5D p ]) in the l.h.s. of eq. (I). Contrary to 
what happened in the Feynman gauge, now such terms do not disappear totally in the 
difference of eqs. (I) and (II). Nor do the terms proportional to f in the r.h.s.'s of Eqs. (I) 
and (II). However, things simplify in the equations for the functions of interest, 5T> a pu k u 
and 5Vl v . 

Consider first the equations for 5T>^ U . When traced over space-time indices, the 
terms depending on 5D^ l \ as well as those coming from the four gluons interactions, 
become identical in eqs. (I) and (II) (at the order of interest). Then, by taking the 
difference of these two equations, we get 



[k ■ D x , 5VT + g j ahc [K8VIX) - (1 - A) (d x SV^k" 

= g(d-l)k v (d^)(d^A)+g(l - r^^KPF-^^A'W, (B.18) 
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which generalizes eq. ( |3.77| ) for arbitrary A. 

Similarly, we shall prove shortly that $D a W satisfies 

(k ■ d x ) (5v; u r) - l g r bc (k ■ a") (8v; u r) + g -f abc k II (A b • 6v c ■ k) 

= ^{KK - g, u k 2 ) {d»A» a ) (d p A) - ^-g 2 r bc (k ■ A b )A^A. (B.19) 

Thus, XdV^k" satisfies the same equation as 5V^ u k u does in Feynman gauge (see eq. ( |3.80| )), 
and, according to eq. (|3.82| ), 

\5V^k»(k,X) = gA;(X)A(k)-k^g a (k,X). (B.20) 



We can use this relation to eliminate the terms dependent on 5V® u k u in eq. ( |B.18|) , which 

k ■ D x , bVl - (1 - A" 1 ) (8Q + g(k ■ A)A'(k))] a 

= g{d-l)k v {d p A v a ){d k p A) + gf abc (k ■ A b )5G c . (B.21) 

This equation shows that 8V> V V — (1 — A -1 ) (5Q + gik ■ A)A'(k)) is independent of A; for 
A = 1, (|B~2TD reduces to eq. flOf ). 

By using the expression QB.20|) in eq. (|3.76|) , we obtain the following expression for 
7 = 7 A - 

J^(k,X) = ^NA(k)A%X) - Nk» (5Vr(k,X) + l -5g\k,X)^ . (B.22) 

It is now straightforward to verify that J g = — (1 — \~ l )J x , where j7" g F is the 
expression obtained for j£ in Feynman gauge, eq. (|3.83| ), while J x is given by eq. QB.14|) . 
Indeed, from eqs. ( [B.14] ) and ( |B.22| ), we see that 

J^ + (l-\- l )J Xa = gNA{k)Al 

-Nk p [5V a v v - (1 - A" 1 ) (5G a + g(k ■ A a )A\k)) + 5Q a } , (B.23) 

where the r.h.s. is, according to eq. ( |B.21| ), independent of A, and equal to (compare, 
in this respect, eqs. ( |B.21|) and ( |3.84| )). Thus, J E + Jn — J g F + j£ , as claimed. 



To complete our proof, we still have to establish eq. (|B.19|) for &T> a k v ' . The technical 
difficulty is the elimination of the terms involving the next-to-leading order correction, 
8D^ ~ g5D^ . In order to do this, we consider the equation satisfied by the derivatives 
of 5D(x,y). First, we obtain an equation for d v y bD ab v by differentiating Eqs. (I) and (II) 
with respect to d y and subtracting them. The l.h.s. of the resulting equation is 

(flg - d 2 y W y 5D; b u ) - (1 - \)d*(d^5Df x ). (B.24) 
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At leading order, 5D^ appears only in the second term: 8^(8% 5 D^) ~ g 2 T 3 D 0l 
i.e., of the order of (8 2 x - 8 2 y ){d v y 8Df v ) « -2{d s ■ d x )(8»5D$ ab ) ~ g 2 T 3 D (recall that 
(d^dgSDp^) = 0). By the same procedure, we get an equation 8^.8 D^, whose l.h.s. is 

{d 2 x - %){%6D*) + (1 - \)dv(d^6D? x ). (B.25) 

Then, eq. ( |B.19| ) for 8T> a lw k v is obtained by subtracting the equations corresponding to 
(P-25|) and ( [B . 24j) . The l.h.s. is, at the relevant order, 

A{d s ■ d x WJD^ ab ) + (1 - \)d%{d^8Df x ). (B.26) 

The second term in eq. ( |B.26| ) is now negligible at this order, as it is ~ g 3 T 3 D . Getting 
the r.h.s. of this equation presents no difficulty A Fourier transform with respect to s 
then yields eq. ( pl9|) for 5V a k v . 



Appendix C 



We prove here that the induced sources rf nA ', eq. (|4.17|) , and j™ d = jt+j^, eqs. ([4.27]) 



and ( |4.39| ), are generated formally by Si n d — S{ + S^, with Sf and Sb given, respectively, 



by eqs. Q5.44| ) and Q5.45| ), if the average fields i/)(X) and A a ^(X) belong to the space 1Z, 



defined by eqs. (|53^)-(|O0|). 



Recall first that, within TZ, the inverse operators 1/v ■ D and l/(v ■ D) 2 which 
appear in eqs. ( |5.44| ) and ( |5.45| ) are unambiguous and represent indifferently the retarded 



or advanced Green's functions, which are equivalent on TZ. For instance, we can easily 
verify, using eqs. (|5.4j ) and (|5.5|) , that 

/ d 4 Y (G ret (X, Y; v) - G adv (X, Y; v)) 4>(Y) = -i \ duU(X, X - vu)^(X - vu) = 0, 

(c.i; 

where the last equality follows from eq. (|5.38| ). This allows us, when doing formal op- 



erations in TZ, to replace 1/v ■ D by the retarded Green's function, which has a simple 
expression in terms of the parallel transporter: 

/ ] r°° 
d 4 Y(X\ -\Y)ip(Y) = duU(X,X -vuU(X -vu). (C.2) 
v ■ D Jo 

This is convenient as it allows us to use nice identities satisfied by the parallel transporters. 
Similarly, we can identify on TZ the operators l/(t> ■ D) 2 and 

G®(X,Y;v) = -Jd 4 ZG ret (X,Z;v)G ret (Z,Y;v) 

= 6(X° - Y°)(X° - Y°) <5 (3) (X-Y - v(X° - Y )) U(X, Y). (C.3) 



59 



Thus, on 7Z, 



d 4 Y(X\- F -\Y)F ia/ (Y)= / duuU(X,X -vu)F^(X -vu)U(X -vu,X).(CA) 

J [y ■ uy Jo 

We conclude that the previous expressions for Sf, eq. (|5.44|) , and for Sb, eq. ( |5.45|) , are 
equivalent respectively to 



S f 



-tuft I ^ I d A X I duip(X)^U(X, X - vu)ip(X - a ) 



(C.5) 



and 



-,S:l^l<rXl dun 



tr {v"F^(X) U(X, X - vu) v v F vX {X - vu) U(X - vu, X)} , (C.6) 

where the trace acts on color indices only. 

Let us prove now that, when the conditions (|5.38|) -( |5.40| ) hold, the induced sources 
are indeed generated from the effective action, according to eqs. ( |5.46|) - ( [5.471) . The first 
equality in eq. (|5.46|) follows immediatly from eqs. (|C.5|) for Sf and ( [4.17D for rf nd . To 
prove the second identity, which relates Sf to j^, we differentiate eq. (|5.44| ) with respect 
to A£(X), using 



6M(X) 



UvD) 



and obtain 



6St 



\Y 2 ) = gv»(Yi 



g^J^v.Jd^Jd^ 



i(v-D) 



\X)t a (X\ 



i(v-D) 



(C.7) 



5A£(X) 



i(v-D] 



\x)i)t a (x\ 



1 



i(vD] 



\Y2MY2 



This is precisely the expression (|5.17|) for (recall that, on 71, G a d v (X, Y; v) = G ret (X, Y; v) 
= (X\l/i(v.D)\Y)). 

We now turn to the verification of eq. ( [5.47] ), which proceeds in three steps. In order 
to do this, it is convenient to use the expression of Sb in terms of the parallel transporter. 

i) We differentiate functionally Sb, eq. ( |C.6| ), with respect to A a ^(X) and write the 
result of this operation as 



5S h = 5S$P + 5Sl 2) . 



(C.9) 
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The first piece comes from differentiating the F^s in eq. ( |C.6| ), 



= l<4f^J d ' X ' duuv»5F; x (X) 
{U ab {X, X - vu) v u F»\X - vu) + U ab (X, X + vu) v v F b vX [X +vu)}. (CIO) 



The second term, 5S^ \ comes from the parallel transporter, 



^Sf ) = 7"IIt L I d ' X I duuv ^ F Hx( x ) $U ab (X, X - vu) v v Fp{X - vu). (C.ll) 



An 



From eq. ([4.14 ), we derive 
5U ab (X,X-vu) 



d Y 



6U ab (X,X-vu) 
5A»(Y) 



5A a JY) 



igu ds U{X,Z)v5A{Z)U(Z,X - vu) 



ab 



(C.12) 



with = T a Al and Z»{s) = X* - u(l - s)v^. We introduce (|CH2D into (|CTTT1) and 
derive, after elementary changes of integration variables, 



5S, 



(2) 



- % 9 7 u p 



f^i [d A x r du r dt in + 1) 

J 4:71 J JO JO 



{v"F^ x (X + vt) [U(X + vt, X) v ■ 5A(X) U(X, X - vu)] v v F b vX [X - vu) j(C13) 



ii) We now derive a new expression for the induced current j£, which will make 
easier the comparaison with SS^/SA^. We note that eq. ( f4.9| ) can be solved in the form 



J Aa {k } X) = gv'k v d p A eff (k) I duU ab (X,X- v'u)F b up (X - v'u) 



(C.14) 



whith ^ = k^/k and A eff (k) = 2N f A(k) + {d- 2)NA(k). The expression ftCTD of the 
induced current j£ is obtained from ( |C.14|) after a direct integration over k. Alternatively, 
one can integrate by parts with respect to k and derive, after some elementary calculations, 



- uj'i / — V,. j du 



Att Jo 



d_ 

dvP 



v^U ab (X, X - vu)F b p {X - vu) 



(C.15) 



Here, v^ = (l,k/ek), but the four components have to be considered as independent 
variables before taking the derivative d/dv p . The integrand in eq. QC.15 ) is further trans- 
formed by writing 



d 

dvP 



U ab (X,X -vu)F b " p (X -vu) 







-u 



+u 



dXP 

U ab (X,Y) 



OXp 



U ab (X,X -vu)F b up (X -vu) 

Fr(X-vu), (C.16) 



Y=X-vu 
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and the second term in the r.h.s. is evaluated by using [22 

D p x U(X,Y)\ y _ x _ vu = igu ^ dss[u{X,Z)v x F xp (Z)U{Z,X - vu) 



(C.17) 



with F xp = T a F xp and Z"(s) = X^ - u(l - s)iA Introducing ([CloD and (|CH7D into 
(|C.15|) and using the identity 



/ 


l - uiv ■ b x ) 


/o 





ab 



U bc {X, X - vu)F^ p (X - vu) = 0, 



we derive the following expression for the induced current: 



with 



#W = £\x)+£\x), 



(C.18) 



(C.19) 



g w (v ■ D x ) - VlM D p ] ab U bc (X, X - vu)v v F v c "{X - vu), 



(C.20) 



and 



3 r dQ r°° r 1 

duu 2 dss 



p 2 J 47T JO Jo 



[/(X, Z) v A F Ap (Z) C/(Z, X - vu) v v F h vp (X - vu). 

J ab 



After using the identity 

U(X,Z)T a U(Z,X) = U ab (Z,X)T b , 

we make a change of integration variables in eq. ( |U.21 ) and finally obtain 

3 /" df2 



(C.21) 



(C.22) 



j^ {X ) = i g tulj^j~duj~dt{u-t) 
[v p F b pX (X - vt) [t)(X - vt, X) T a U(X, X - 



be 



,Ff{X-vu)}. (C.23) 



iii) It is now straightforward to verify that, for gauge fields in TZ, the following 
equalities are true 



(C.24) 



6SP/6A£(X) = j« a (X), SSP/SA&X) = j^ a (X). 
The first equality results from eqs. ( |C10| ) and ( |C.20| ), by recalling that 

5F^(X) = b^A v {X)-D u 5A^X), (C.25) 
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and by using the identity 

J duuU ab (X, X - vu)F* u (X - vu) = J duuU ab (X, X + vu)F* u (X + vu), (C.26) 



which follows from the condition ( 5.4U| ). As for the second equality ( fj.24 ), we note that 
( CT) and (|CTD imply 

, . Q /■ ^/O /-oo /-oo 

6S®/6A£(X)-jV>*(X) = -ig-^lj ^ / + 

{w p F p 6 A (X + t;t) [f/(X + t»t,X)T a f7(X,X-m)] ;ic t» i ,F ( r A (X-t;M)}. (C.27) 

Again, this vanishes because of the conditions (|5.39|) and ( |5.40|) . This completes our proof. 
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FIGURE CAPTIONS 



Fig.l: 

Diagrammatic illlustration of the equation of motion (2.18) for the average gauge field 
(An(x)). The thick line represents the gluon propagator in the presence of the background 
field (Ap), i.e. with all possible mean field insertions included. The cross represents the 
external current j)x(x), while the circled cross represents the induced current j™ d (x) whose 
detailed contributions in terms of 2-point and 3-point functions are displayed in (b) (see 
eq.(2.22)). The contributions of the ghost mean fields are not drawn. Quark, gluon and 
ghost propagators are represented respectively by straight, wavy and dashed lines. 

Fig.2: 

Diagrammatic illustration of eqs. (2.19) for the fermionic field (ip{x)) (a), and its con- 
jugate {4>(x)) (b). The thick lines represent quark propagators in the presence of the 
background gauge field (A^). The crosses stand for the external sources r\ and fj, and the 
induced sources are represented explicitly by the appropriate 2-point functions (see eq. 
(2.23)). 

Fig.3: 

Lowest order contribution to the off equilibrium fermion propagator (a). The blob rep- 
resents a gauge field insertion and carry a four momentum P ~ gT. Typical radiative 
corrections are shown in (b). The contributions in (b) are negligible compared to (a) 
when k ~ T, but are of the same order as (a) when k ~ gT. 

Fig. 4: 

Illustration of the approximate equations of motion for K^(x,y), eq. (3.22) (a), and eq. 
(3.23) (b). Thick lines represent propagators in the presence of the background gauge 
field (Afj) (see Fig. 5a for the quark propagator and Fig.6a for the gluon one). 

Fig.5: 

Illustration of the approximate equation of motion (3.24) for S(x, y). (a): quark propaga- 
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tor with gauge field insertions; (b): corrections induced by the fermionic fields; the thick 
lines represent propagators in the presence of the background gauge field (A^). 

Fig.6: 

Illustration of the approximate equation of motion (3.27) for D(x, y). (a):gluon propagator 
with gauge field insertions; (b): corrections induced by the fermionic fields (only the 
contribution due to the K function is shown); the thick lines represent propagators in the 
presence of the background gauge field (A^). 

Fig.7: 

The induced current j md calculated from the 2-point functions obtained in section 3.3. In 
the first two lines are represented the various contributions to j A = jf+jg+jc+Jn- m the 
last two lines are represented, respectively, the contributions jf and to = jf + jf, 
as obtained by joining the end points x and y in Figs. 5 and 6. The ghost mean field 
contributions are not drawn. 

Fig.8: 

The induced fermionic source rf nd calculated from the 2-point functions obtained in section 
3.3, that is, by joining the end points x and y in the diagrams in Fig. 4. 



67 



